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Abstract 

We study the partition function ZQ(^nk,k){QT''^) of the Q-state Potts model on the 
family of (non-planar) generalized Petersen graphs G{nk,k). We study its zeros in 
the plane {Q,v) for 1 < A; < 7. We also consider two specializations of ZQ(^nk,k)^ 
namely the chromatic polynomial PG{nk,k)iQ) (corresponding to v = —1), and the flow 
polynomial ^G{nk,k){Q) (corresponding to v = —Q). In these two cases, we study 
their zeros in the complex Q-plane for 1 < k < 7 . We pay special attention to the 
accumulation loci of the corresponding zeros when n — )• oo. We observe that the 
Berker-Kadanoff phase that is present in two-dimensional Potts models, also exists 
for non-planar recursive graphs. Their qualitative features are the same; but the main 
difference is that the role played by the Beraha numbers for planar graphs is now played 
by the non- negative integers for non-planar graphs. At these integer values of Q, there 
are massive eigenvalue cancellations, in the same way as the eigenvalue cancellations 
that happen at the Beraha numbers for planar graphs. 



Key Words: Potts model; non-planar graphs; Beraha conjecture; generalized Petersen 
graphs; transfer matrix; Berker-Kadanoff phase. 



1 Introduction 



The two-dimensional (2D) Q-state Potts model [41,58] is one of the most studied models 
in Statistical Mechanics. Despite many efforts over more than 60 years, its exact free energy 
and phase diagram are still unknown. The ferromagnetic regime of the Potts model is the 
best understood case: exact (albeit not always rigorous) results have been obtained for the 
ferromagnetic-paramagnetic phase transition temperature T(.{Q) for several regular lattices, 
the order of the transition (continuous for < Q < 4, and first order for Q > 4), the 
phase diagram, and the characterization in terms of conformal field theory (CFT) of the 
corresponding universality classes. (See e.g., Ref. [3].) 

The antiferromagnetic (AF) regime is less understood. This is partly because, in contrast 
with the ferromagnetic regime, universality cannot be expected to hold in general. Investiga- 
tions must therefore proceed on a case-by-casc basis. For instance, the free energy is known 
exactly along some curves of the phase diagram (Q, T) (where T is the temperature), for cer- 
tain regular 2D lattices [3,4]. One of these curves belongs to the ferromagnetic regime, and 
it can be identified with the ferromagnetic-paramagnetic phase-transition curve. It might be 
tempting to infer the very existence of a phase transition from this (partial) solubility of the 
model along that curve. One well-known example of the invalidity of such an inference is 
provided by the zero-temperature limit of the triangular-lattice Q-state Potts antiferromag- 
net [5,6]. Although its free energy is exactly known for all values of Q e R,^ the system is 
known to be critical only in the interval Q e [0, 4], and disordered for Q e (— cxo, 0) U (4, oo). 

In three dimensions (3D) there are no known exact results for the Q-statc Potts model. 
Most numerical results come from series expansions and Monte Carlo simulations: see e.g., 
Refs. [23,37, and references therein]. In the ferromagnetic regime, we expect a critical curve, 
which is second order for Q ~ 2, and first-order for Q — 3. There was an important 
controversy in the late 80 's and early 90 's about the precise nature of the Q = 3 transition 
because of its relation with QCD: the four-dimensional SU(3) lattice gauge theory should be 
in the same universality class as the 3D ferromagnetic 3-state Potts model [37]. 

The Q-state Potts model at temperature T can be defined on any (undirected) finite 
graph G — (V, E) with vertex set V and edge set E. On each vertex i e y , we place a 
spin that can take Q distinct values: Uj e {1,2, ... ,Q}. These spins interact through the 
Hamiltonian [41] 



(ij)eE 

where 5ij is the usual Kronecker delta, and J is a real coupling constant that is proportional 
to 1/T. The partition function is defined as usual as: 



Notice that initially, Q is a positive integer Q>2, and J is a real number. The ferromagnetic 
(resp. antiferromagnetic) regime corresponds to J > (resp. J < 0). 

^The Q-state Potts model can be defined for non- integer values of Q using the Fortuin-Kasteleyn repre- 
sentation explained below [c.f. (1.3)]. 




(1.1) 




(1.2) 
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Fortuin and Kasteleyn [22] have shown that the partition function (1.2) can be rewritten 

as 

Zg{Q.v) = 5] ^l^'lg'=(^') , (1.3) 

E'CE 

where the sum runs over the 2l^l subsets E' C E, with k{E') being the number of connected 
components (including isolated vertices) in the spanning subgraph {V,E'), and v is the 
temperature-like parameter 

V = e-^-1. (1.4) 

We now can promote (1.3) to the definition of the model, which permits us to consider the 
parameters Q and v as arbitrary complex numbers, because (1.3) is a polynomial in both. 
The Fortuin-Kasteleyn (FK) representation of the Q-state Potts model (1.3) is equivalent 
to the Tutte polynomial studied by graph theorists [55] after a change of variables. In 
terms of v, the ferromagnetic (resp. antiferromagnetic) regime corresponds to v > (resp. 
— l<v< 0). The unphysical regime corresponds to v < —1 (i.e., complex J), and the 
behavior of the Potts model in this regime is less well understood than that of the AF 
regime. 

There are several interesting particular cases of Zo[Q,v). One first example is the chro- 
matic polynomial 

Pg{Q) = Zg{Q,-1), (1.5) 

which corresponds to the zero-temperature limit in the AF regime J — )■ —oo. The restriction 
of this polynomial to Q a positive integer has indeed an interpretation as a coloring problem 
[16]: PaiQ) gives the number of proper Q-colorings of the graph G. A proper Q-coloring of 
a graph G = {V, E) is a coloring of the vertices in V such that for each edge e = {hj) G E 
its endpoints i.j are not colored alike. 

The second example is the flow polynomial 

MQ) = (-i)i^ig-i^i^G(g,-Q). (1.6) 

Again, the restriction of this polynomial to integer Q has a combinatorial interpretation: it 
gives the number of nowhere Zq-Aows on G. If P is an additive Abelian group of order Q, 
a P-flow on G is a function (p: E F, so that each edge e ^ E is associated to a variable 
0(e), subject to the constraint that these variables are conserved at each vertex i E V, given 
an arbitrary orientation of the edges in E. A no-where zero P-flow is a P-flow such that 
0(e) 7^ for all edges e E E [36,56,60]. If P is a finite Abehan group of order Q, then 
the number of nowhere zero P-flows depends only on Q (not on the specific structure of the 
group P), and it is in fact the restriction to Q G N of a polynomial in Q called the flow 
polynomial $g(Q) [54]. We choose P = Zq as our group of order Q. 

These two polynomials are related for planar graphs. This relation is based on the duality 
transformation for the Q-state Potts model on a planar graph G [59] : 

Zg{Q, v) = gi^i-i^i-i ^1^1 Zg* {Q,v*), (1.7) 

where G* is the dual graph of G, and v* is the dual of v 

vv* = Q. (1.8) 
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Then for v — —1, v* — —Q, so that (1.7) reduces to 

Pg*{Q) = Q^aiQ), (1.9) 

The duahty relation (1.9) can also be understood combinatorially: the flow along an edge 
e & E of G determines the color difference between the two faces (vertices of G*) that are 
adjacent to e. The conservation of the Z^-fiow at each vertex i ^ V ensures that the color 
differences thus defined add up to zero upon encircling i. A definite Q-coloring of G* is 
obtained from these color differences upon fixing the color of one reference vertex; this is 
responsible for the factor of Q appearing in (1.9). 

For non-planar graphs there is no known relation whatsoever between $g(Q) ^-nd Pq* (Q)- 
The evaluation of Zg{Q,v) for a general graph is a hard problem, viz., at least as de- 
manding as the determination of the chromatic or the flow polynomials. In particular, the 
determination of the coefficients of the chromatic or fiow polynomials of a general graph (in- 
cluding bipartite planar graphs) is #-P hard [38, Proposition 2.1], and the same thus holds 
in general for the computation of the coefficients of the full partition function Zg{Q,v). For 
recursive families of graphs, the partition function Zg„{Q,v) of any member of the family 
Gn (composed by n identical layers of size m each) can be computed from a transfer matrix 
T and certain boundary condition vectors u and v: 

ZgAQ, v) = u{mY ■ T(m)" • v{m) , (1.10) 

where t denotes the transpose. The computation time thus grows as a polynomial in n; but it 
docs however still grow exponentially in m. Indeed, there arc similar formulas for computing 
Pg and $g directly from the corresponding transfer matrices. The partition function (1.10) 
can be written as a sum over the eigenvalues Xj of the transfer matrix T with appropriate 
amplitudes aj: 

dimT 

i=i 

where both the amplitudes and eigenvalues are algebraic functions of both Q and v. 

From (1.10) one can compute the partition function Zg„ of the Q-state Potts model (or 
any of its specializations Pg„ or $g„) on any member Gn of a recursive family of graphs. 
In practice, this family of graphs will be a 2D strip graph of some regular lattice of width 
m and length n with some boundary conditions.^ In this case, we will simplify the notation 
by writing Z^^y^np, where the subscripts a,P = F,P denote free and periodic boundary 
conditions in the respective lattice directions. For 3D "slab" graphs, we use ^maxp^xn^ for 
the partition function of a slab of section m x p and thickness n. As the partition function 
is a polynomial in Q and v, we can study its roots, e.g., by fixing Q (resp. v) to a physical 
value and then obtaining the roots of the resulting one-variable polynomial in the complex 
v-plane (resp. Q-plane). This programme has been done for 2D strip graphs of the square 
and triangular lattices with free and cylindrical boundary conditions [17,18]. Indeed, there 

^We adopt Shrock's [51] terminology for boundary conditions of 2D strip graphs: free (rrip x np), cylin- 
drical (mp X np), cyclic (mp x np), and toroidal (mp x np). Here the first dimension (m) corresponds to the 
transverse ("short" or space-like) direction, while the second dimension (n) corresponds to the longitudinal 
("long" or time- like) direction. 
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are several ways to study the zeros of Zq^'- for instance, considering the plane {CliV) where 
both variables take real values. This approach will produce the "phase diagram" of the 
model for finite transverse size. 

The above study is hard to carry out in practice. In fact, most studies are concerned 
with the zeros in the complex Q-plane of the chromatic polynomial Pq (or chromatic zeros 
of G). In a series of papers [27-29,31,46,47, and references therein], we have studied 2D 
strip graphs of the square and triangular lattices with free, cylindrical, cyclic, and toroidal 
boundary conditions. For 3D lattices, the authors of Ref. [45] computed the chromatic roots 
of "slabs" with some small transverse sections (with free and periodic boundary conditions) 
and several values of the longitudinal size n (with free boundary conditions). Finally, the 
roots of the flow polynomial (or flow roots) have been studied for some 2D strip graphs and 
boundary conditions in Ref. [20]. 

From these zeros it is very hard to obtain infinite- volume quantities (m, n — )■ oo), as they 
have strong finite-size-scaling (FSS) corrections. Moreover, one must face the well-known 
difficulty that the limits lim^=„_).oo (Fisher limit) and lim^_^oo hnin^oo (van Hove limit) may 
give different results in some parts of the parameter space. In this paper we consider the 
latter limit, which has two major advantages: 1) the first limit, lim„_^oo, is convenient (and 
accessible in polynomial time) within the transfer matrix approach, and 2) the resulting FSS 
dependence in the second limit, lim^_j.oo, is much weaker than when considering the m and 
n limits simultaneously. In other words, we first compute, for fixed (and small) values of m, 
accumulation sets of partition function zeros in the infinite- length limit n — > oo, and study 
the infinite-width limit m — )■ oo subsequently. 

According to the Beraha-Kahane-Weiss (BKW) theorem [10-13,53], these zeros accumu- 
late along certain limiting curves (when the two eigenvalues that are largest in modulus 
become equimodular) , and around isolated limiting points (when there is a unique dominant 
eigenvalue and its amplitude vanishes) . By computing the eigenvalues \j and amplitudes aj 
of the transfer matrix T (1.11), we are able to obtain the exact values of the relevant physical 
quantities in the limit n — )■ oo. Their FSS corrections are found to be smaller than when 
both m and n are finite. The extrapolation to the true infinite- volume limit m ^ oo, using 
standard FSS techniques, therefore becomes more precise by employing this order of limits. 

When looking at the real chromatic zeros for 2D strip graphs of the square and triangular 
lattices with free, cylindrical and cyclic boundary conditions [27,28,31,46,47, and references 
therein], we observe (as many authors in the literature did in the past!) that there exist real 
accumulation points of these real chromatic roots. ^ This observation dates back to Beraha [9] 
who made a conjecture about the possible values of Q that can be accumulation points for 
real chomatic roots. As explained by Saleur [49] there are two statements of the Beraha 
conjecture that do not seem to be fully equivalent: The first one is contained in Beraha's 
PhD Thesis [9] (the statement is taken from Ref. [11]): 

Conjecture 1.1 (Beraha vl) Among the limits of all (not necessarily recursive) families 
of chromatic polynomials are found the numbers of the form: 

= 2 + 2COS— = 4cos2-, (1.12) 
n n 

•^The same is true for the full partition function zeros of strip graphs of the square and triangular lattices 
with free and cylindrical boundary conditions for fixed physical values of t; > —1 [17, 18]. 
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for any positive integer n. Notice that Bi = lim Bn — 4. 

This conjecture has a shghtly different form in Baxter's book [3] (the statement is taken 
from Ref. [6]): 

Conjecture 1.2 (Beraha v2) Some of the real zeros of chromatic polynomials of planar 
graphs should, in the limit of the graph becoming large, occur at points in the sequence Bn 
(1.12) with integer n>2. 

Finally, Jackson [25] states the Beraha conjecture in another slightly different way: 

Conjecture 1.3 (Beraha v3) There exists a plane triangulation with a real chromatic root 
in {Bn — e,Bn + e) for all n>2 and all e > 0. 

Remeirk. In the context outlined above, our understanding is that the Beraha conjecture 
is a statement about the accumulation points of chromatic roots in the infinite-size limit 
of planar regular graphs. Conjecture 1.2 claims that generically the Beraha numbers are 
limiting points for the chromatic zeros of families of planar graphs. It is clear that integer 
Beraha numbers B2, B3, B^, Bq can actually be chromatic zeros of planar graphs (e.g., K4 
has all these numbers as chromatic roots: Pk4,{Q) = Q{Q — 1)(Q ~ 2)((5 — 3)). However, 
non-integer Beraha numbers (except perhaps Bio) cannot be chromatic zeros of any planar 
graph [46, Corollary 2.4]. Furthermore, Bio is not a chromatic zero of any plane near- 
triangulation [46, Proposition 2.3(c)], nor of any plane triangulation. 

Let us now consider a strip graph of a (not necessarily planar) regular lattice with periodic 
boundary conditions along the longitudinal direction. To treat this case within the transfer 
matrix formalism, we have to keep track of the bottom- and top-row connectivities. As 
explained in detail in Refs. [28-30] , the transfer matrix only acts on the top-row connectivity, 
so the full transfer matrix has a block-diagonal form, each block corresponding to a different 
bottom-row connectivity. Indeed, there are blocks of the top-row state connected to blocks of 
the bottom-row state; each of these structures will be called a link (or bridge). As the transfer 
matrix cannot increase the number of links of a given connectivity state, if we order the state 
appropriately, the transfer matrix takes an upper-triangular form. Therefore, the eigenvalues 
can be obtained from those diagonal blocks corresponding to a given number i of links. The 
corresponding amplitudes are non-trivial and can be inferred either from a combinatorial 
reasoning applied to the upper-triangular decomposition of the transfer matrix [19,43,44], 
from quantum field theory [42,50], or from representation theoretical considerations [24]. We 
now have three cases: 

(1) The strip is planar with free boundary conditions along the transverse direction. This 
is the geometry of an annulus. In this case, the links cannot interchange their positions, 
so the group acting on them is the trivial group E consisting only of the identity. The 
transfer matrix commutes with the generators of the quantum algebra Uqsl(2), with 
VQ = q+q^^- In this case the amplitudes ae are given by Chebyshev polynomials of the 
second kind [19,40,50]. The numerical evidence [28] shows that there are accumulation 
points at the Beraha numbers (1.12).^ This is compatible with the CFT predictions of 
the pattern of eigenvalue dominance [50]. 

"^We find [27,31,46,47] that the Beraha numbers are accumulation points also for strip graphs of the 
square and triangular lattices with free and cylindrical boundary conditions. 
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(2) The strip is planar, but periodic boundary conditions are imposed along the transverse 
direction. The resulting geometry is that of the torus, which is obviously non-planar. 
In this case, the links can be cyclically interchanged across the strip. Therefore, the 
group acting in this situation is the cyclic group C^. The representation theory of this 
group leads to very different expressions for the amplitudes [44, Equations (1.2)/(1.3)] 
(see also [42]) involving number theoretical functions. The fact that the I links can now 
wind around the periodic transverse direction with some momentum p implies that the 
amplitudes a^^p depend on both parameters. In this case the numerical evidence [29], 
and CFT arguments for the pattern of eigenvalue dominance, shows that there are 
accumulation points of real zeros at (^ = 0, 1,2, 3, 4. 

(3) The strip graph itself is non-planar. In this case the links can be interchanged in 
any way, so the group describing this situation is the full symmetric group S^. The 
amplitudes a^^x obtained from the representation theory of this group are certain poly- 
nomials that depend on I and on the irreducible representation A of S(, [24] . The best 
numerical evidence comes from the study of the flow polynomial of the generalized 
Petersen graphs [30]. In this case there are accumulation points of real roots at non- 
negative integer values oi Q = 0,1,2,.... The present paper extends this evidence to 
higher values of Q, and to the Potts model partition function in general. 

Remark. Obvious a planar strip graph is a special case on a non-planar one; yet the results 
(3) do not apply to the cases (1) or (2). This is because the transfer matrix enjoys more 
symmetries (i.e., has a larger commutant) in the planar case, and this must be taken into 
account in the corresponding representation theory. In the same vein, it might be that some 
special classes of non-planar strip graphs lead to higher symmetries than generic non-planar 
strips (see e.g. [21]) and hence to modifications of the results (3). In this paper we shall 
give compelling evidence that the generalized Petersen graphs provide an example of generic 
non-planar strip graphs. 

In practice, one does not see all the real accumulation points expected from the above 
scenario. In fact, only those within the Berker-Kadanoff (BK) phase [49,50] can be observed. 
This is a massless phase with algebraic decay of correlations throughout which the temper- 
ature parameter v is irrelevant in the renormalization group (RG) sense. In the generic 
case, one expects the following picture [33] of the phase diagram of the 2D Potts model on 
a given regular lattice in the {Q,v) plane. We start with the ferromagnetic critical curve 
vfm{Q) > 0; the transition is second order for < Q < 4, and first-order for Q > 4 [2]. 
Along this curve, the thermal operator is relevant, and becomes marginal in the limit Q ^ 
of spanning trees. The analytic continuation of vpu into the AF regime is a critical curve 
denoted VBKiQ) < with < Q < 4. Along Vbk the thermal operator is irrelevant. We can 
think of this curve as a renormalization-group attractor, whose basin of attraction is the BK 
phase. This phase is bounded by two curves v±{Q). The upper one v+ is usually identified 
with the AF critical curve. ^ 

This picture has been verified [28, 29] for the square and triangular lattice. For the 
square lattice we know the exact form of these curves i'fm(Q) = +VQj 'VBKiQ) — —VQ, 

^By definition this is the phase transition curve in the v < region that is closest to the infinite- 
temperature limit V = 0. 
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and v±{Q) — —2 ± — Q [2-4]. Exactly v — v±{Q) we find a different type of critical 
behavior [32]. For the triangular lattice we know that v-pM, vbk-i and are the tree branches 
of the equation [3, 7]: 

v^ + Zv"^ = Q. (1.13) 

If we denote as Q2(— 1) the first point where v+ crosses the line v = —1, then one can only 
see the accumulation points of chromatic zeros in the interval [0,(52(— !)]• For the square 
lattice, we know that Q2{—1) = 3, while for the triangular lattice, Q2{—1) ~ 3.81967 [5,6].^ 
It is important to stress that for 2D Potts models the BK phase does not exist when Q is one 
of the Beraha numbers i?„ (1-12). This is due to the fact that at these values of Q, some of 
the amplitudes vanish, and there are eigenvalue cancellations, which give rise to a different 
physics. These phenomena were already found in the study of the chromatic polynomial 
of strip graphs of the square and triangular lattices with cyclic [28, Section 6.4] boundary 
conditions. 

If we move to lower values of v < — 1 (i.e., in the unphysical phase), the BK phase extends 
to a w-dependcnt value Q2{v). For any given lattice, we denote by Qmax the maximum value 
of Q2{v) that one can attain within the BK phase; obviously then Qmax > Q2(— !)• One 
has Qmax = 4 for both the square and the triangular lattices [2,4,7,28,29]. Moreover, the 
representation theory of Uqsl{2) implies [40] that one has Qmax < 4 for any 2D Potts model. 

Remcirk. The main features of the BK phase in 2D, including the special role of the Beraha 
numbers, have recently been confirmed by the independent technique of graph polynomials 
[34,35]. 

To attain larger values of Qmax, one must thus study families of non-planar graphs [like in 
point (3) above]. In this paper, and motivated by our findings in Ref. [30], we will consider 
the generalized Petersen graphs G{m,k) (sec Section 2). These are generically non-planar 
graphs that can be built in a recursive way. We expect the generalized Petersen graphs to 
be representative for generic non-planar strip graphs. Moreover, they are devised to have 
large values of Qmax, so that we gain access to a vast swath of the BK phase. 

We find for this family of graphs that: 

• The qualitative picture of the phase diagram in the {Q, v) plane for the (non-planar) 
generalized Petersen graphs agrees well with that of 2D Potts models. In particular 
we find a BK phase. 

• The value of Qmax is numerically found to be large, Qmax ^ 12.4(1). 

• The set of non- negative integers {0, 1,2,...} plays the same role as the Beraha numbers 
did for 2D Potts models. 

• At these integer values of Q, we find amplitude vanishing and eigenvalue cancellations. 
This leads to a scenario that ressembles that of Uqsl{2) symmetric 2D Potts models. We 
give details of the inclusion/exclusion compensation of eigenvalues for Q = 0,1, 2, 3, 4. 

^However, in Ref. [31] a re-analysis of Baxter eigenvalues given in Refs. [5,6] yielded the value Q2{—1) = 
Bx2 « 3.73205. 
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• In the phase diagram in the {Q,v) plane we find regions where there is not a single 
dominant eigenvalue, but a pair of complex conjugate eigenvalues. According to the 
BKW theorem, we expect that these regions will contain a dense set of partition- 
function zeros. In previous studies all zeros were found to accumulate on points or 
curves, so this is a genuinely new feature. 

The plan of this paper is as follows. In Section 2 we describe the most important prop- 
erties of the generalized Petersen graphs G{m,k). In Section 3 we summarize the basic 
properties of the transfer- matrix formalism applied to the family of graphs G{nk,k), and 
in Section 4 we will show our results for the partition-function zeros and limiting curves in 
the {Q,v) plane. In Section 5, we will analyze in detail what happens to the eigenvalues 
and amplitudes when Q is a non-negative integer. In Sections 6 and 7 we will consider the 
zeros and limiting curves in the complex Q-plane of two specializations of Zc{Q, v), namely 
the fiow and chormatic polynomials, respectively. Finally in Section 8 we discuss our results 
and make some general conjectures about the role of Beraha (resp. non-negative integers) 
for planar (resp. non-planar) graphs. In Appendix A we will show another example of non- 
planar graphs (the 3D simple-cubic graph of section 2x2) with similar properties to those 
of the graphs G{m, k). 



2 The generalized Petersen graphs 

Wc shall consider the family of graphs G{m, k) called generalized Petersen graphs and 
defined as follows: Let m, k be positive integers such that m > k. Then G{m, k) is a cubic 
graph (i.e., all vertices have degree 3) with 2m vertices denoted ip and jp for p = 1, 2, . . . , m: 
i.e., 

V{G{m,k)) = (2.1) 

The edge set consists of 3m edges {ipjp), {ipip+i), {jpjp+k), for p = 1,2,..., m, and with all 
indices considered modulo m: i.e., 

E{G{m, k)) = {(ipJp), (ipip+i), (jpjp+k) | 1 < p < m} . (2.2) 

Note that G{m, k) is simple for m ^ 2k; but it has double edges when m = 2k. These graphs 
were introduced by Watkins [57]. As an example, G(12,4) can be drawn as follows: 
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The graphs G{m, k) are non-planar for all pairs (m, k), except for the case (3, 2), and the two 
sub-families {p, 1) and (2p, 2) with p > 1. More properties of these graphs can be found in 
Ref. [30, and references therein]. In this paper, we will focus on the sub-family of generalized 
Petersen graphs G{nk, k) with n>2. Transfer-matrix methods can be efficiently applied for 
this family (see next section). 

3 Potts model transfer matrix 

We wish to evaluate Zc{Q, v) (1.2) / (1.3) for G = G{nk, k) by a transfer matrix construc- 
tion. Contrary to an often repeated but false statement, evaluating Zg{Q,v) by a transfer 
matrix construction is possible for any graph G, and does not require G to consist of a num- 
ber of identical layers [8]. However, when G does have a layered structure — as is the case 
here — Za{Q, v) can be computed by the repeated application of the same transfer matrix. 

The construction of the transfer matrix for the partition function of the generalized 
Petersen graphs G{nk,k) was done in detail in Ref. [30]. So here we will just summarize 
the main results. The first step is to redraw the graph G{nk, k) in such a way that it is 
composed by n identical layers, each one containing 2k vertices, and with periodic boundary 
conditions in the longitudinal direction. Each layer has width L — k + 1. (See [30, Figure 1].) 
All edges now link vertices within the same layer, or in two adjacent layers. 

We will work on a space of partitions of the set {0, 1, 2, . . . , A;, 0', 1', 2', . • • , k'}, where 
{0, 1, 2, . . . , A:} (resp. {0', 1', 2', . • • , k'}) belong to the top (resp. bottom) row of the strip. 
Given a partition, a link is a block that contains vertices belonging to both the top and 
bottom rows. The number of links in a given partition will be denoted by £. 

The transfer matrix Tl adds one layer on top of the strip. It can be written in terms of 
vertical Vj and horizontal Hjj operators: 

Hij = I + v^ij (3.1a) 
Vi = vI+Di (3.1b) 

where Jj^ is a "join" operator that amalgamates the blocks containing the points i.j, and Dj 
is a "detach" operator that removes point i from its block and turn it into a singleton (if i 
was already a singleton it provides a factor Q to that term). Then the full transfer matrix 
can be written as 

Tl = Hoi (n^oHo'O (n^M = HoiVoHo2VoHo3...VoHofeVfeVfc_i...Vo, (3.2) 

\i=k / \i=0 ) 

where the rightmost operators act first. 

As the transfer matrix only acts on the top-row connectivity state, it is independent of the 
bottom-row state. So it can decomposed into a diagonal block form, each block corresponding 
to a different bottom-row state. All these blocks give the same eigenvalues, so we can choose 
a simple one. As the transfer matrix cannot increase the number of links ^, then if we order 
the states in an appropriate way, the transfer matrix can be written in an upper-triangular 
block form. Indeed, the eigenvalues can be obtained from its diagonal blocks each of 
them corresponding to a given value of the number of links 0<£<L = A;-|-1. Finally, given 
a block with i links, these links can be interchanged in any possible way, so the relevant 
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group is the symmetric group S^. Then it turns out [24] that the relevant eigenvalues come 
from the irreducible representations A of of dimension dim A. If wc choose as our basis the 
appropriate linear combination of states, then the transfer matrix T/^^^ takes a diagonal block 
form, where each block T ^^^^x corresponds to a distinct irreducible representation A e S^. 
Then, the partition function for G{nk, k) can then be written as a sum of ordinary traces: 

fc+i 

ZG{nk,k)iQ,v) = J] J]autr(Ti,^,;,)". (3.3) 

e=o xeSe 

The amplitudes ae^x are polynomials in Q given in Ref [24, Proposition 3.24]: 



dim A ^ ^ 



oie,x 



i=0 



Yl{Q-i-X,_,), (3.4) 



where the representation A G 5*^ is considered through its corresponding Young diagram, 
F(A) = (Ai, A2, . . . , A^), where Aj is the number of boxes in the i'th row. If there is less than 
i rows in ^(A) the expression is of course padded with zeros. 

We have symbolically computed all relevant blocks T^^x forl<A;<6, 0<£<A; + 1, and 
A G Se. Therefore, we could in principle compute the partition function for all generalized 
Petersen graphs G{nk, k) with 1 < k < 6 and arbitrary n > 2. It is worth stressing that 
Eq. (3.3) holds true irrespective of whether some eigenvalues happen to be identical or 
not. However, it will turn out useful in practice to have a "complete' decomposition of the 
partition function, in the sense that there are no equalities among eigenvalues (except for 
what is accounted for in the multiplicities a^ x)- 

As for the particular case of the flow polynomial $G(nfc,fc) [30], we have found that for 
1 < A; < 6, the blocks T^x with L = k + 1 have an upper-triangular block form: 



k+i,e,\ — 




''k+i,e,x 

-p(nt) 

' k+i,e,x 



(3.5) 



where D^+i^^ ;^ is a diagonal matrix with all its diagonal terms are equal to the "trivial" eigen- 
value /ifc,fc+i — v'^^, and diagonal block containing the non-trivial eigenvalues 
(hence the superscript "(n^)"). Notice that for £ = 0, there are no trivial eigenvalues: i.e., 
dim Dk+ifi = 0; while for £ = A; + 1, all eigenvalues arc trivial: i.e., dimT^"|*| ^^-|^ a ~ 0- "^^^ 
dimension of the diagonal matrix D^+i^e^x is the same as for the flow-polynomial case [30]. 

Finally, we have also observed that for 1 < /c < 6, the block T^*^} ^ ^ contains only 2A; -|- 1 
non-trivial eigenvalues iik,k,s of multiplicity dim A for all irreducible representations X & Sk- 
These eigenvalues come from a reduced transfer matrix T^*^^^ ^. 

Therefore, for 1 < A; < 6 the partition function for the generalized Petersen graphs 
G{nk, k) (3.3) can be re-written as 

fc-i 

ZGink,k){Q, ^) = E E tr(Tgl_,^,)- + /3, tr(Tg^)'^ + ^,^,v'-^ , (3.6) 

where the coefficient is given by the sum 

A = E «^,AdimA, (3.7) 
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and the coefficients 7^+1 are given by 

k 

Ik+i = h+i + X] ^^'^ ^™ ^k+i,t,\ ■ (3.8) 
i=\ XeSe 

The expressions for the needed and 7^+1 can be read from [30, Eqs. (3.28)/(B.ll)]. 

We have proven (3.6) by exphcit computation for 1 < /c < 6. We conjecture that (3.6) is 
also true for each k >7. This conjecture indeed holds true in the particular case of the flow 
polynomial [30]. 

If we denote by i^k,e,\s (resp. i^k,k,s) the eigenvalues of T^"*|^^_;^ (resp. T^'^*^), then (3.6) 
can be written as 

fe-l Nt,{e,X) 2k+l 

ZGink,k){Q^^) = I]5Z"^,A Yl l^kAX,s + l3kJ2l^lk,s + ^k+lv'''', (3-9) 
e=0 XGSe s=l s=l 

where Nk{i,X) is the number of non-trivial eigenvalues corresponding to i links {0 < £ < 
A; + 1) and the irreducible representation A e Se- The expression (3.9) can be regarded as a 
"complete" decomposition of the partition function in the sense that for a fixed value of k, 
all the eigenvalues appearing in (3.9) (namely, fJ'k,e,x,s, f^k,k,s, and v^*^) are all distinct. This 
fact was proven for 1 < /c < 6 by explicit computation of the eigenvalues. Furthermore, we 
conjecture that (3.9) holds true with all eigenvalues being distinct also for k >7. The total 
number of distinct eigenvalues is 6, 20, 113, 755, 5568, 43975 for /c = 1, . . . , 6. 

As each eigenvalue in (3.9) is characterized by the number of links i, we can use this 
number as a sector label. Therefore we will use the expression "sector i" (or sector of £ 
links) as a shorthand for the set of non-trivial eigenvalues with number of links equal to £ 
{0 < £ < k). The sector £ = k + 1 will correspond to the trivial eigenvalue iik,k+i = v^''- 

4 Numerical results 

In this section we provide information about the practical procedure we have followed to 
compute the relevant blocks of the transfer matrix. We then show our results concerning the 
phase diagram for each value of k. Finally, by comparing these finite-A; results, we attempt 
to gain some insight about the thermodynamic limit k ^ 00. 

4.1 Practical procedure 

We have written a MATHEMATICA script to compute the symbolic transfer matrix Tl 
(with L = k + l and for 1 < A; < 6) using ideas similar to those already explained in [28-30] . 

The first goal is to obtain the relevant diagonal blocks Tk+i,e,x of the transfer matrix 
Tl- We first fix the number of finks £ (0 < ^ < A; + 1) and a bottom- row state compatible 
with the chosen value of £. We then determine the basis of connectivity states. The result 
indeed docs not depend on the chosen bottom-row state. We now choose an irreducible 
representation A G S*^ of dimension dim A. The relevant diagonal block T^+i^^^a is obtained 
by taking as our basis vectors those linear combinations of the "standard basis" with the 
appropriate symmetries under S£. To decrease the CPU time needed for the computation. 
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we extract the trivial eigenvalues Hk^k+i — vhy looking for a upper-block diagonal form like 
in (3.5). 

Once the basic blocks t[,"|*| ^ ^ are obtained, we compute the traces tr(T^'^| ^ y)"^ using 
again MATHEMATICA. We then reconstruct the partition function using (3.6)-(3.8). 

We have checked the partition functions for small values of n by computing the partition 
function using maple. Indeed, from the full partition function ZG(nk,k){Q,'v) we can obtain 
its flow-polynomial specialization ^G{nk,k){Q)- These polynomials were also compared to 
those obtained by other methods in Ref. [30]. We find a perfect agreement in all cases. 

We have symbolically computed all relevant blocks T^+i/^x ior 1 < k < 6, < £ < k + 1, 
and \ E Si- We have also found the "complete" decomposition of the partition function for 
these values of k. However, the actual computation of the traces is quite demanding: for 
A; < 5 we were able to compute the partition function for values n < 10. For A; = 6, we could 
only compute n = 1, 2. The partition functions for n = 3, 4 were obtained using MAPLE. For 
k = 7 we could not obtain the transfer-matrix blocks; but we managed to get the partition 
functions up to n = 4 using again maple. 

4.2 Results for the generalized Petersen graphs 

Our goal is to study the "phase diagram" of the Q-state Potts model on the generalized 
Petersen graphs G{nk, k) in the {Q, v) plane. From a physical point of view, this is the most 
natural approach to the problem. We will focus on the antiferromagnetic/unphysical region 

V < 0, as we are interested in studying a possible BK phase in this model. 

To study the partition-function zeros in this plane, we took several lines: Q = —pv, 

V = —pQ (with p = 1,2,3), and v = — 1, — 2, . . . , — 6. The line v = —Q (resp. v = — 1) 
corresponds to the fiow-polynomial (resp. chromatic-polynomial) subspace. Along these 
lines, the two- variable partition function ZG{;nk,k){Q-iV) reduces to a single- variable polyno- 
mial with integer coefficients. We then find its zeros by using the program MPSOLVE by 
Bini and Fiorentino [14,15]. These lines are depicted as dot-dashed lines in Figures 1 and 2, 
and the zeros are also represented as circles and squares in the same figures. Although 
they exhibit some systematic features, these zeros are generally too scattered for any firm 
conclusions to be drawn. 

A much better approach is to compute the limiting curves Bfc when n — >■ oo in the (Q, v) 
plane. We have used the direct-search method [46] to compute Bk- But in some areas of 
the plane, this method did not work properly. The reason was that there are regions where 
the dominant eigenvalue is actually a pair of complex conjugate eigenvalues belonging to the 
same sector. These regions are labeled by an asterisk in Figures 1 and 2. 

For A; = 7, 8, 9, 11, even though we could not obtain the relevant transfer-matrix blocks, 
we were able to numerically compute (parts of) the corresponding limiting curves by using 
a code written in C. For 2 < A; < 6 the result of this program coincides with the previous 
symbolic computation. 

By looking at Figures 1 and 2, we see that for each fixed value of k (with 1 < A; < 7) 
there are two distinct regions in the corresponding phase diagram: 

(a) This region is the closest to Q = and ends approximately at Q ~ A; -|- 1. It contains 
vertical lines at even integer values of Q, and has isolated limiting points at odd 
integer values of Q. It is bounded above and below by simple curves, and the dominant 
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Real zeros of Z^^^^ ^^(Q.v) Real zeros of Z^^g^^ gj(Q,v) 




Figure 1: Real zeros of the Potts-model partition function and limiting curves Bk in the 
plane {Q,v) for the generalized Petersen graphs G{nk,k) with k = 1 (a), k = 2 (b), k = 3 
(c), and k = 4 (d). For each panel, we show the real zeros along several lines (depicted 
as doted-dashes lines) of the type Q = —pv {p = 1,2,3), v = —pQ {p = 1,2,3), and 
V = —1, —2, —3, —4, —5, —6. For each value of k, the zeros correspond to the generalized 
Petersen graphs G{9k,k) (black □) and G{10k,k) (red o). Each region is labeled with the 
sector the dominant eigenvalue belongs to (e.g., i = 1). The asterisk in the sector label 
(e.g., i = 1*) means that there is a pair of complex conjugate dominant eigenvalues in that 
region. The label i = (1^) in panel (c) means that the dominant eigenvalue comes from the 
completely anti-symmetric irreducible representation (1, 1) = (1^) of 5*2. 
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Figure 2: Real zeros of the Potts- model partition function and limiting curves Bk in the plane 
{Q,v) for the generalized Petersen graphs G{nk,k) with k = 5 (a), k = 6 (b), and k = 7 
(c). In (a) the zeros correspond to the generalized Petersen graphs G{45,k) (black □) and 
G(50, k) (red o). For each value of k in (b) and (c), the zeros correspond to the generalized 
Petersen graphs G{3k,k) (black □) and G{Ak,k) (red o). The rest of the notation is as in 
Figure 1. For k = 7 we could not obtain the more involved structure for large Q ^ 8.2. 
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eigenvalues are simple ones, hence real valued. Close to Q ^ 0, the dominant eigenvalue 
belongs to the sector i — 1, and every time we cross one of these vertical lines the value 
of i increases by one unit. 

(b) This region starts approximately at Q ^ k + 1, and it has a more involved structure. 
For k = 1 there is a single outward branch starting at (Q, v) = (4, —2) and extending to 
infinity. For 2 < A; < 6, we find two boundary curves that extend outwards to infinity. 
Inside this region there are sub-regions where the dominant eigenvalue is actually a 
pair of complex conjugate eigenvalues (in most cases they belong to the sector ^ = 0). 
For A; = 3, we even find some regions where the dominant eigenvalue comes from 
the completely anti-symmetric representation (1,1) = (1^) of the i — 2 sector (see 
Figure 1(c)). 

Above [resp. below] both regions the dominant eigenvalue comes from the i — [resp. 
£ = (A; + 1) mod 2] sector. 

Region (a) looks qualitatively similar to the BK phase we have found for the 2D Potts 
models on the square and triangular lattices. Inside the region enclosed by the limiting curve 
Bk, we find that for each integer p > up to some maximum value Pmax, there is a "phase" 
corresponding to an eigenvalue belonging to the fully symmetric representation A = {£) of 

with i = p + 1 links. This phase contains all values of Q in the interval Q G (2p, 2p + 2). 
Therefore, there are phase transitions at Q = 2p for p > in an interval {v_{p),v^{p)), 
in which eigenvalues from the £ = p and £ = p + 1 become equimodular. As k increases, 
the maximum value Pmax increases, so we find more phases in the BK phase. Even though 
Q2(— 1) < 3 is rather small (see Section 7), as in the 2D case, we now find vertical lines at 
Q- values as large as Q = 8 for k = 7 {p = 4), which is the maximum value of k for which we 
have been able to produce the full "regular part" — region (a) — of the phase diagram. 

We emphasize that the presence of vertical lines of limiting points extending from (Q) 
to v+{Q) provides clear evidence for the extent of the BK phase. In particular, the perfect 
verticality means that the temperature parameter v is RG irrelevant for V-{Q) < v < Vj^{Q), 
and the physics is thus controlled by an RG fixed point in the interior of the BK phase. Ideally 
one would like to corroborate the identification — and characterization — of the BK phase by 
verifying the algebraic decay of correlation functions inside the BK phase. In the 2D case 
this is a rather simple matter, since standard CFT results relate the strip free energies to 
critical exponents. In Section 4.3 we pursue a more modest goal, namely to give numerical 
evidence that the correlation length ^ is proportional to k inside the BK phase and tends to 
a constant outside it. Since k roles as the effective width in the strip geometry underlying 
the transfer matrix formalism, this implies that the BK phase is characterized by ^ oo in 
the thermodynamical limit As — > oo, and so is critical indeed. 

In the phase characterized by £ links, corresponding to the interval Q e {2{£ — l),2£), the 
amplitude corresponding to the fully symmetric representation A = {£) is [c.f., (3.4)] [24,30] 

atlt) = ^(Q-2^+l)n(Q-i). (4.1) 

Therefore, it vanishes at the mid-point Q — 2£ — \ oi the above interval; hence, this is an 
isolated limiting point. 
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Figure 3: Limiting curves Bk in the plane {Q,v) for the generahzed Petersen graphs G{nk, k) 
with k = 1 (black), k = 2 (gray), k = 3 (light blue), k = 4 (green), k = 5 (orange), k = Q 
(pink), k = 7 (violet), k = 8 (red), k = 9 (navy blue), and A; = 11 (black). We label some 
regions with the sector to which the dominant eigenvalue (for large enough k) belongs (e.g., 
i = 3). The horizontal v = —1 (resp. inclined v = —Q) dot-dashed line corresponds to the 
chromatic (resp. flow) polynomial subspace. 



Remark. Notice that the families G{n, 1) and G{2n, 2) are planar. Indeed, in Figures 1(a)- 
(b), we see that there are isolated limiting points at the non- integer value Q = Br, in 
agreement with the Beraha conjecture. This can be explained by the fact that the amplitude 
for the i = 2 sector is not the generic one, but 72 for = 1, and (32 for k = 2. (cf., [30, 
Eq. (3.28)]). Indeed, both amplitudes are equal to /32 = 72 = — 3Q + 1, and P2{B^) = 0, 
as expected. 

Figure 3 shows all the limiting curves put together. If we look at the bottom part of 
these curves {Q = 0,v < —4), the value of k increases as we move outwards. We notice the 
following empirical observations: 

• As k increases, the BK phase enlarges and more i sectors are visible in the unphysical 
region. 

• The upper boundary curve converges quickly to a limit. In particular, its slope at 
Q = seems to be equal to —1/3.'' Its crossing with the horizontal line v = —1 

^For each value of fc = 1, . . . , 6, we have fitted the data closest to the origin to a polynomial Ansatz. For 
odd values of k, we obtained in all cases estimates compatible with a slope equal to —1/3 (e.g., —0.333333(1) 
for k = 1,3, and —0.3333(1) for k = 5). For even values of k, we obtain estimates that seem to converge 
to -1/3 from below (e.g., -0.347841(1), -0.339861(1), and -0.3368(1) for k = 2,4,6). A power-law fit to 
these three data points gives a slope of « —0.329, not far from the odd-fc estimate. Therefore, the numerical 
data suggest that the slope of the curve v+{Q) at Q = is exactly —1/3, as claimed. 
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(depicted as dot-dashed line in Figure 3) defines the value of Q2{—1). (See Section 7.) 
We can identify this curve with the AF critical curve v+{Q) for this model. 

• The lower boundary curve V-{Q) converges slowly to some limit. The inclined dot- 
dashed line in Figure 3 corresponds to the flow polynomial. The convergence is not as 
fast as for (52 (—1)- (See Section 6). 

• Using the rough estimates for the hmiting curves for A; = 9, 11, we see that Qmax ^9.4, 
which is far larger than the maximum value Qmax = 4 obtainable for 2D Potts models 
with either cyclic and toroidal boundary conditions. We can obtain a more accurate 
estimate of Qmax by looking at the line v = —4 (see Figures 1-2). This line seems to 
cross the "regular" part of the BK phase at its more distant point from the origin, 
and has the nice property that the eigenvalue that is dominant on the high-Q side 
of the crossing comes always from the £ = sector for all values of k we can access 
numerically: i.e., 3 < A; < 10.^ These results are displayed in Table 1. As a prevention 
against parity effects, we have fitted the odd- and even-Zc data subsets separately. We 
have used several Ansatze: i.e., power law and polynomials of different degrees in 1/k. 
From the dispersion of the estimates for Qmax, we conclude that our preferred estimate 
is Qmax 12.4(1), as the limiting value of the crossings along the line v — —A only 
provides a lower bound for the true Qmax- 



k 


Qcik) 


3 


3.5918146982 


4 


5.5837328676 


5 


6.6418204973 


6 


7.7533996189 


7 


8.3495427730 


8 


8.8740090544 


9 


9.2660844176 


10 


9.5702423022 



Table 1: Values of Qc{k) along the line v — —A for the generahzed Petersen graphs G{nk, k) 
for 3 < A; < 10. 



4.3 Criticality of the Berker— Kadanoff phase 



The BK phase in 2D Potts models is a critical phase: it corresponds to a massless phase 
with algebraic decaying of correlations. In this section we want to study whether the BK 
phase for the non-planar Petersen graphs G{nk, k) is critical or not. One way to achieve this 
goal is to consider the inverse correlation length (or mass gap) given by the ratios 



^-\Q,v) = log 



(4.2) 



^For fc = 11, Arnoldi's method did not converge, probably due to a complicated eigenvalue nearby. 
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^ along Q=3/2 for Petersen 



^ along Q=3/2 for Petersen 




-8 -6 -4 -2 

V 

(a) 

^ along Q=7/2 for Petersen 



8 -6 -4 -2 

V 

(b) 

^2 ^ along Q=7/2 for Petersen 




Figure 4: Inverse correlation lengths ^^^{Q,v) for Q = 3/2 (a,b), and Q = 7/2 (c,d). The 
index i = 1 (resp. i = 2) corresponds to the the ratio of the first (resp. second) sub-dominant 
eigenvalue and the dominant one [c.f., (4.2)]. We show the inverse correlation lengths 
{i = 1,2) for the Petersen graphs G{nk,k) with k = 1 (black), k = 2 (gray), k = 3 (light 
blue), k = 4 (green), k = 5 (orange), k = 6 (pink), k = 7 (violet), k = 8 (red), and k = 9 
(navy blue). 
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^ along Q=3/2 for Petersen ^ along Q=3/2 for Petersen 




-0.6 -0.4 -0.2 -0.6 -0.4 -0.2 



V V 

(a) (b) 

Figure 5: Inverse correlation lengths ^^^{Q,v) for Q = 3/2 and z = 1 (a) and i = 2 (h) for 
the Petersen graphs G{nk, k). The color code is as in Figure 4. We show in both panels the 
curves for 1 < /c < 9. 



where /x* is the dominant eigenvalue, and jji with i = 1,2,... is the z-th sub dominant 
eigenvalue (i.e., > > |/i2| > • • •)• Figure 4 we have shown the values of for 
i = 1,2 versus the parameter v at fixed values of Q = 3/2, 7/2.^ Note that for Q = 3/2 we 
cross the region of the BK phase dominated by the i = 1 sector, while for Q = 7/2, we cross 
the nearby region dominated by the i = 2 sector. We expect to find a disordered phase for 

V G {v+{Q),0], the BK phase for v G {v-{Q),V-^-{Q)), and another phase for v < V-{Q). The 
critical/non-critical nature of the last two phases needs to be determined. 

Let us now start with the first region v G (f+(Q),0]. Figure 4 does not give much 
information about this regime; so we provide a zoom of this region in Figure 5 for (5 = 3/2. 
Now the behavior of the different curves as k increases is clear: they tend to some finite limit 
> 0. This limit curve is essentially attained for /c > 3 and v > —0.3. We conclude 
that in this region there is a finite gap in the thermodynamic limit, and therefore, this is 
a non-critical phase. For Q = 5/2 the plots are very similar. To summarize, the region 

V G (f+((5),0] is a non-critical (disordered) phase. 

Let us now consider the last region v G {—oo,v^{Q)). The most relevant part of this 
regime is displayed in Figure 6 for Q = 3/2. The plot of ^ displayed in Figure 6(a) show 
clearly that this inverse correlation length tends to zero as k increases. Therefore, one could 
naively conclude that this also a massless phase. There are two almost identical eigenvalues, 
one coming from the i = sector, and the other from the i = 1 sector. But the plot of 
in Figure 6(b) show that the second correlation length is finite. For k > A and v > 7.2 
all curves collapse almost perfectly. Therefore, we expect that as k tends to infinity, these 
curves will tend to some finite limit ^2^^(f) > 0. The plots for Q = 5/2 are very similar. 
Therefore we conclude that the region v G (— oo,f_((5)) is also a non-critical phase. The 



^We have also made computations for Q = 11/2, but the results are qualitatively very similar to those 
for Q — 7/2, so we refrain from displaying the corresponding plots. 
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^ along Q=3/2 for Petersen 



^2 ^ along Q=3/2 for Petersen 




Figure 6: Inverse correlation lengths ^~^{Q,v) for Q = 3/2 and i = 1 (a) and i = 2 (h) for 
the Petersen graphs G{nk, k). The color code is as in Figure 4. We show in both panels the 
curves for 1 < /c < 9. 



two most dominant eigenvalues for finite values of k will tend to a common function, so the 
mass gap in this phase is given in the thermodynamic limit by ^2^^. 

Finally, we focus on the BK phase v G {v-{Q),V-^.{Q)). Figure 4 shows that in this region 
the curves are rather flat around some non-zero value for k > 3, and as k increases, this 
value becomes smaller. This is an indicator that as k increases the mass gap tends to zero in 
this phase. Therefore, we expect that this is a critical (massless) phase (as for the 2D Potts 
model). In 2D models, CFT predicts that the mass gap for a critical theory defined on a strip 
graph of width L and cylindrical boundary conditions behaves like = Ai/L + o(L~^), 
where Ai is a universal constant, proportional to the critical exponent that governs the 
algebraic decay of an appropriate correlation function. For the non-planar recursive strip 
graphs studied here, there is no such strong link to CFT, but it is nevertheless true that 
k roles as the effective width of the strip in the transfer matrix formalism [30]. Critical 
behavior will therefore result if the mass gap for the non-planar Petersen graphs G{nk, k) 
behaves asymptotically like 1/k. To determine whether this is the case, we have plotted in 
Figure 7 the quantity k^~^ versus v for fixed values of Q = 3/2, 7/2. 

Even though there are obvious finite-size effects, for Q = 3/2 it looks like the curves tend 
to some limiting curve from below. However, for Q = 7/2,^^ this convergence is from above. 
This means that there is a change of sign in the dominant finite-size-scaling correction term. 
The best convergence is clearly attained for Q = 7 /2 and i = 1. We have tried to fit the 
data for Q = 3/2 and 7/2 and i = 1, 2 to an Ansatz 

k^r\Q,v) = g.{Q.v) + ^^^^, (4.3) 
where the second term intends to mimic the finite-size-scaling corrections. For the four 



For Q — 11/2, the finite-fc curves seem to tend to some limiting curve from above. 
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Figure 7: k times the inverse correlation length, k^l^{Q, v), with (5 = 3/2 (a,b) and Q = 1 /2 
(c,d) for the Petersen graphs G{nk,k) in the region v G {v-{Q),v+{Q)) for z = 1 (a,c) and 
i = 2 (b,d). For Q = 3/2, we show in panels (a,b) the curves for 1 < A; < 9. For Q = 7/2, 
we only show in panels (c,d) those curves with 3 < < 9, as for smaller values of k the line 
(5 = 7/2 crosses the limiting curves Bk at their non- regular part. The color code is as in 
Figure 4. 
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cases we have considered here, we find reasonable fits leading to smooth functions gi{Q, v) 
as functions of v. Therefore, we conjecture that in the BK phase 

Cr\Q,v) - ^^ + o(0- (4.4) 

This implies that the BK phase is critical in the sense that it is massless. This coincides 
with the properties of the BK phase found for 2D Potts models. 

In Section 4.2, we have shown that the limiting curves Bk contained vertical lines for 
even integer values of Q = 0,2, . . . and v in the interval v G (f -(Q), 'y+(<5))- The perfect 
verticality of these lines implies that the temperature is an irrelevant operator in the RG 
sense. Moreover, in Section 4.3, we have shown that inside the region v e {v-{Q),V-^{Q)) 
for fixed values of Q, the mass gap vanishes in the infinite- /c limit. This implies that this 
region is critical. Therefore, both conditions that define the BK phase (as we know from 2D 
Potts models) are met also for the Q-state Potts model on the Petersen graphs G{nk, k). 

We conclude that the BK phase exists also at least for recursive families of non-planar 
graphs, and a role analogous to the one played by the Beraha numbers for planar graphs is 
now played by the non-negative integers for non-planar graphs. The natural question is to 
know, in the non-planar case, whether there arc also cancellations among eigenvalues and 
vanishing of amplitudes when Q is a non-negative integer. This question will be answered 
in the next section. 



5 What happens for integer values of Q7 

In this section we will consider the partition function for the generalized Petersen graphs 
ZQ{nk,k){Q-i'^) when Q takes a non-negative integer value. 

Let us fix A; > 3. For each < £ < /c, let us denote by £i the set of non-trivial eigenvalues 
f^k,e,\,s associated to the sector of i links. Moreover, let S^+i be the one-element set containing 
only the trivial eigenvalue Hk,k+i = ^"^^ ■ In the same way, let £(,^x denote the set of non-trivial 
eigenvalues iik,i,x,s corresponding to the irreducible representation A of the group S^. 

When (5 is a non-negative integer, some of the non-trivial eigenvalues of T^*^^^ ;^ (for 
< ^ < A;-l) or T^^^l^ J, (for £ — A;) may reduce to the trivial eigenvalue Hk,k+i — 

These 

eigenvalues are excluded by definition from the sets and E^^x defined above for < ^ < /c. 

Finally, let us denote by Xk,i,\ the net non-trivial contribution to the partition function 
Za{nk,k) [c-f- (3.9)] of the sector with i links and the irreducible representation A of Sf. 

Xk,e,x = Yl l^h,\s- (5-1) 

s 

The prime on the sum in the r.h.s. of (5.1) means that we exclude, for a particular non- 
negative integer value of both 1) those non-trivial eigenvalues of T^*^^ ^ ^ that become 
identical to the trivial one v'^^, and 2) those eigenvalues in S^^x that are exactly cancelled 
by the contribution of other sectors. For £ = k, wc use the eigenvalues Hk,k of the reduced 
matrix T^!^| ^ instead, and A plays no role. For £ — 0,1, the sub-index A is also superfiuous 
and can be dropped. 
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We can obtain the eigenvalues of each transfer matrix Tk+i,e,\ by computing its char- 
acteristic polynomial symbolically using MATEMATICA. This procedure can be carried out 
for k = 2,3,4; but for k = 5, this is not feasible due to the large dimension of some of the 
transfer matrices. We have observed that the characteristic polynomials obtained in this way 
can not be factorized over the integers. (This is not true in general when Q takes integer 
values, as shown below.) Indeed, if we find that two transfer matrices have characteristic 
polynomials with a common factor, then we can deduce that the eigenvalues coming from 
the common factor appear in both matrices. On the contrary, if the characteristic polynomi- 
als of two transfer matrices have no common factor, then they do not contain any common 
eigenvalue. For k — 5 we have numerically computed the eigenvalues to high precision (i.e., 
100 digits) using mathematica. For these cases, we could check the conjectures obtained 
from the smaller values for k. 

We summarize our findings for the smallest integer values by using the exact results for 
A; = 3,4, 5 (which correspond to non- planar generalized Petersen graphs). 

5.1 Q = 

This case is trivial, as Zg{0,v) = for any graph with at least one vertex. Therefore, 
there should be a massive cancellation of eigenvalues and amplitudes, so that the net sum is 
always zero. 

First of all, it is clear from (3.4), that only the fully anti-symmetric representations of 
Si give a non-zero contribution. In fact, if we represent by (1^), the fully anti-symmetric 
representation of S^, then 

«£,(iO(0) = (-1)'- (5.2) 
The two non-generic amplitudes are also [30]: 

Pk(0) = 7fc(0) = (-l)^ (5.3) 

We find that Sq G Si. i.e., the set of i = eigenvalues is a proper subset of those of the 
next sector. For the higher sectors (except the last one), we find that Se \ Sg-i C Sg+i. In 
words, the non-trivial eigenvalues of the sector with i hnks that do not belong to the sector 
with £ — 1 links, are a proper subset of the set of non-trivial eigenvalues of the sector with 
£ + 1 links. The last step corresponds to £ = k, and it reads Sk \ £k-i = ^- (Recall that 
Si for < i < k contains non-trivial eigenvalues, while £k+i contains the trivial eigenvalue 
t;2^) 

When Q — 0, some of the eigenvalues iJ,k,e,x in (3.9) reduce to the trivial one ijLk,k+i — v'^'^- 
These new trivial eigenvalues appear with some multiplicity rik/ > in sectors 1 < £ < k 
(resp. 2 < £ < k) for odd (resp. even) values of k (e.g., when k = 4, 714 2 = 'n.4,3 = 3, and 
n4^4 = 1). But their net contribution always cancels the last term on the r.h.s. of Eq. (3.9): 

k 

J2nkA-^Y + lk+i{0) = 0. (5.4) 

e=i 

The above inclusion-exclusion patterns, in addition to the values taken by the amplitudes, 
imply that a// eigenvalues cancel out exactly, giving the expected result ZG{nk,k){^,v) = 0. 
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5.2 Q = l 

This case is also trivial, as for any graph G = {V,E), Za{l,v) = (1 + f)''^'. This comes 
from the FK representation (1.3) of the Q-state Potts model. Recall that for the generalized 
Petersen graph G{nk, k), \E\ = 3nk. 

The generic amplitudes Oie,x{l) are non-zero only for the following cases: 

• i — with CKo = 1. 

• i = 2 and the representation (2) with 0:2,(2) = — 1- 

• For each £ > 3, the representation (2, with q;^,(2,i^-2) = (—1)^"^. 
The non-generic amphtudes [30] take the values 

= {k-l){-l)'-\ 7fc+i = (-!)'• (5.5) 

The pattern for the non-trivial eigenvalues is as follows: For i — 0, we have the decom- 
position Sq — {/io}U£^0' with /io = (1 +?;)'^'^. Then an inclusion- exclusion pattern similar to 
the (5 = case occurs: E'l C £"2, S-t \ Se-i C Se+i for 3 < £ < A; — 1, and Sk \ £k-i = 0- 

There is a slight difference with the Q = case: the non-zero amplitude for i = k is 
{k — 1)(— 1)'''^^, so that the corresponding (non-trivial) eigenvalues in Ek contribute to the 
partition function with this amplitude. But these eigenvalues also appear in the preceding 
sector i = k — 1 with a degeneracy exactly equal to /c — 1; therefore, these two contributions 
cancel out in the partition function. 

The trivial eigenvalue v'^'^ cancels out in all cases 3 < A; < 5. Unfortunately, we have 
been unable to find a simple pattern for the multiplicities and the amplitudes in this case. 

Therefore, all eigenvalues cancel out exactly, except the simple eigenvalue — {1 + v)^'' 
coming from the £ = sector. This eigenvalue gives the right value of the partition function 
when Q — 1. Using (5.1), we can write the partition function as 

ZG(nk,k){^,v) = Xk,o = (1 + ^;)'"^ (5.6) 

5.3 Q = 2 

This case corresponds to the Ising model on the family of graphs G{nk, k). The generic 
amplitudes q;^,a(2) (3.4) are non-zero only for the following cases: 

• £ = with tto = 1. 

• I — 1 with a\ — 1. 

• 1 — 2 and the representation (2) with q;2,(2) = — 1- 

• £ = 3 and the representation (3) with Q!3,(3) = — 1. 

• For each £ > 4 and the representations (3, 1^~^) and (2^,1^~^), with 0:^,(3,1^-3) = 
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The non-generic amplitudes take the values [30] : 

13k = (fc'-3fc + l)(-l)^ 7fc+i = -1- (5.7) 

The pattern for the non-trivial eigenvalues is more involved than for the previous cases. 
In the sector £ = 0, we find two types of non-trivial eigenvalues: Sq = 8'q\J 8q. The first 
sub-set contains \6q\ — 2^ eigenvalues, and these eigenvalues do contribute to the partition 
function ZG(nk,k)i'^,v)- The second subset is contained in {£q C £3); but the amplitudes 
of these two sectors are equal in absolute value with distinct signs, so they give a zero net 
contribution to ZG{nk,k)i.'2j v). In the sector £ = 1, a. similar phenomenon occurs: £1 = £[U£'^, 
with = 2*^ contributing eigenvalues, and the rest cancel out with part of the i — 2 sector: 

C £2- The non-trivial eigenvalues in £2 \ £'{ and £^3 \ £q also appear in the higher sectors 
f > 4, and their net contribution is always zero. We also find for A; = 4, 5 and sector £ = k, 
2k + 1 distinct eigenvalues with amplitude The same 2k + 1 non-trivial eigenvalues also 
appear at lower sectors with multiplicities that make the their net contribution equal to zero. 

The trivial eigenvalue v'^'^ cancels out in all cases 3 < k < 5. But its pattern does not 
look simple, and probably has parity effects. 

Therefore, all eigenvalues cancel out exactly, except the 2'^"'"^ eigenvalues in the set £qU£'i. 
This is expected, as on each layer of G{nk, k) there are exactly k + 1 vertices, and hence the 
corresponding transfer matrix in the spin representation has dimension 2*^"'"^. Therefore, the 
partition function reads: 

ZG{nk,k)('^,v) = Xk,0 + Xk,l, (5.8) 

where each Xk,e contains 2^^ non-zero eigenvalues. This formula looks like the expression 
found for the RSOS representation of the 2D Ising model [26] . 

5.4 Q = 3 

The situation is now more involved, but the result is rather simple: 

ZG{nk,k){3, V) = Xk,0 + 2Xk,l + Xfc,2,(l2) , (5.9) 

where there are [3^/2] 3^^, and [3^/2] non-zero eigenvalues in the above "characters", 
respectively. In total, there are 2x3*^ contributing eigenvalues in the set £q[J£[[J£2, where we 
have spht each eigenvalue set £k — £'^[j£'j^, as in the case Q — 2. The final formula resembles 
the formula found for the RSOS representation of the 3-state Potts model in 2D [26]. 

We also find for k = 4,5 the following patterns for the other eigenvalues (whose net 
contribution to ZG{nk,k){'^,'v) cancels out): in the sector £ = 1, £i\£[ C. ^3,(3)- In the sector 
£ — 2, the only non-vanishing amplitude is 0:2,(12) = 1, while for sector £ — 3, there are two 
non-zero amplitudes q;3^(3) = —2 and ck3^{2,i) — —1- We then find that £^2 \ ^2 ^3,(2,i)- Let 
us note that Q = 3 is inside the phase characterized hy £ — 2 links; but from this sector 
most eigenvalues cancel out. 

The trivial eigenvalue v^^ cancels out in all cases 3 < A; < 5. But its pattern does not 
look simple, and probably has again parity effects. 

5.5 Q = 4 

The analysis of the eigenvalues for Q = 4 is more complicated (in fact, for /c = 3 we do 
not find any cancellation at all!). For A; = 4, 5 we do find cancellations; in particular, all 
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eigenvalues associated to negative generic amplitudes a(,^\ cancel out, so that we conjecture 
the following formula: 

ZGink,k){4:, V) = Xk,0 + 3Xfc,l + 2Xfe,2,(2) + 3Xfc,2,(l2) + Xfe,3,(l3) 

even] , (5.10) 

where I[A] is the indicator function for the event A {I[A] = 1 if A is true, and zero, otherwise). 
We cannot determine the value of the constant A (or even if it is really a constant), as we 
only have one data point {k — 5). 

The non-trivial eigenvalues associated to £ = 3 and A = (3) (with 03,(3) = ~2) cancel out 
with those appearing in ^^2,(2)- i-c, ^3,(3) C ^^2,(2)- For k = 5, wc find 11 distinct eigenvalues 
for £ = k. These eigenvalues appear also for {£, A) = (1, (1)), (4, (3, 1)), (4, (2, 1^)); and their 
net contribution is zero. The other 305 eigenvalues that appear in ^^4,(3,1) (resp. £^4,(2,12)) 
cancel with those appearing in ^2,(12) (resp. £^3^(13)). In summary, all eigenvalues associated 
to negative amplitudes cancel out (except the trivial one for even values of k), and there are 
only five non-trivial contributions to the partition function. (Indeed, this is in sharp contrast 
with the infinite sum obtained using the RSOS representation for the 4-state Potts model 
on a planar graph [26]). 



5.6 Conjectured behavior for integer Q 



The above results for Q = 0,1,2,3,4 suggest the following conjecture for the partition 
function of the generalized Petersen graphs ZG(nk,k){Q, v) (3.6) when Q takes a non-negative 
integer value. 

Conjecture 5.1 Let us fix Q to a non-negative integer value N. Then the partition function 
of the generalized Petersen graphs (3.6) reduces for any k > N to: 



N-2 



ZG{nk,k)(N,v) = ^ 



E 



+ Xfc,JV-i,(i)^-i +Pfe(-^)^^ 



2nk 



(5.11) 



where we have defined the "characters" Xk,e,x in (5-1), Pk is the number of the trivial eigen- 
value contributes to the partition function, and the second sum is over all irreducible repre- 
sentations X E Se with a positive value of the amplitude ai^\ at Q = N . 



Remcirks. 1. Notice that o;jv-i,(i)Jv-i = 1. 

2. The function Pk{N) is simple for small values of < A?" < 4: 



jO ifO<iV<3 , 

b l[k even] if A/ = 4 



3. For any k > N, we find the following sum rule: 



AT-l 



J]«,,A7Vfe(£,A) + pfe(7V) = Ar*^+i, (5.13) 
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where {i, A) is the number of contributing non-trivial eigenvalues to the partition func- 
tion ZG(nk,k){N,v) and coming from the sector with i links and the representation A G S£. 
Note that (5.13) is compatible with the fact that, for Q integer, the transfer matrix can 
alternatively be represented in terms of Potts spins. In this representation the partition 
function with periodic longitudinal boundary conditions is an ordinary matrix trace, and no 
eigenvalue cancellations can take place. Therefore, the number of eigenvalues (counted with 
multiphcities) is Q'^'^^ indeed. 



6 Flow polynomial for the generalized Petersen graphs 

In Ref. [30] we focused on the particular case of the flow polynomial of the generalized 
Petersen graphs G{nk, k) for 1 < k < 7. Using the eigenvalues found there, our present 
goal is to compute the accumulation sets of flow zeros in the limit n — > oo. In particular 
we consider the limiting curves Bk for 1 < < 7. We have computed these curves for 
k = 1,2, ... ,7 using the direct-search method, as described in [46]. The results are shown 
in Figures 8 and 9 in the complex Q-plane. 

We have noticed several empirical patterns in Figures 8 and 9. The first one is related 
to the number of outward branches of the limiting curve Bk, and the dominant eigenvalue 
in each asymptotic sector 

Conjecture 6.1 Fix k > 1. Then, the limiting curve Bk has exactly 2k outward branches 
extending to Q — oo, with asymptotic angles arg Q — O^j where 

Ok,j = (j-^)p forj = l,...,2k. (6.1) 
Moreover, the dominant eigenvalue is Hk+ifi,* in the asymptotic regions 

a ^ ^ a t / J = 1, 3, . . . , 2/c - 1 if k is even 

9k,j < argQ < 9k,j+i for < f u ■ ^^'^^ 

\ J — 2,4, ... ,2k if k IS odd 

while in the other asymptotic sectors the dominant eigenvalue is iJ>k+i,i,{i),*- 

We may also conjecture that, as A; — )■ oo, the limiting curves Bk (without the outward 
branches) converge to some curve Boo- This is what Figure 10 suggests. Indeed, outside this 
curve, the outward branches seem to get denser as k increases. Therefore, we conjecture 
that the set of accumulation points of the flow-polynomial zeros for the generalized Petersen 
graphs G{nk, k) is dense in the whole complex Q-plane, except in the interior of the curve 
Boo- If this conjecture is true, it implies in particular that flow roots of non-planar cubic 
graphs are dense in the complex plane, except possibly in some bounded region near the 
origin. Note that a similar result holds (by duality) for the flow roots of planar graphs [53], 
but without the restriction to cubic graphs. 



^"'^This is exactly the same empirical behavior found for the limiting curves Bm associated to the chromatic 
zeros of the family Sm,n when n — >■ oo [48, Conjeture 7.1]. The graph Sm,n can be regarded as a square-lattice 
grid with m columns, n rows, free boundary conditions in the longitudinal direction, and special boundary 
conditions in the transverse direction: wc introduce two extra vertices such that every vertex on the leftmost 
column of the grid is connected to one of them, and every vertex on the rightmost column is connected to 
the other one. 
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Complex zeros ^q^^ 
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Complex zeros *(j(4ii4) 




Figure 8: Complex zeros of the flow polynomial and limiting curves Bk in the complex 
Q-plane for the generalized Petersen graphs G{nk,k) with = 1 (a), /c = 2 (b), /c = 3 
(c), and k = A (d). For each value of k, the zeros correspond to the generalized Petersen 
graphs G{10k,k) (black □) and G{20k,k) (red o). The labels (e.g., i = 1) show the sector 
the dominant eigenvalue belongs to. We only label the regions that have a non-empty 
intersection with the real axis. 
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Complex zeros *Q(5n 5) 



Complex zeros *Q(6n 6) 
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Figure 9: Complex zeros of the flow polynomial and limiting curves in the complex Q- 
plane for the generalized Petersen graphs G{nk, k) with = 5 (a), /c = 6 (b), and k = 7 (c) 
The black squares (□) correspond to the zeros of G{50, 5) (a), G(60, 6) (b), and G{63, 7) (c) 
The red circles (o) correspond to the zeros of G(100, 5) (a), (7(120, 6) (b), and G(119, 7) (c) 
Labels are as in Figure 8. 
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Limiting curves for *(j(iikk) 




Re(Q) 

Figure 10: Limiting curves Bk for the flow polynomial in the complex Q-plane for the gen- 
eralized Petersen graphs G{nk, k) with k = 1 (black), k = 2 (gray), k = 3 (light blue), k = 4 
(green), k = 5 (orange), k = 6 (pink), and k = 7 (violet). Labels are as in Figure 8. 



If we restrict to the cases with k > A (see Figures 8(d) and 9), we notice some regularities 
about the isolated limiting points and the dominant eigenvalues in the non-asymptotic region 
of the complex Q-plane. If we denote by Qc*'' the largest real value of Q where Bk crosses 
the real Q-axis, then we conjecture that: 

Conjecture 6.2 Fix k > A. Then, the dominant eigenvalue in the regions that intersect the 
real Q-axis is 



(6.3) 



Aife+i,i,{i),* for Re Q G (-oo, 2] 
yUfc+i,2,{2),* for ReQ G [2,4] 
/i* = /^fc+i,3,{3),* forReQ e['i,Qf\k)] 

A*fe+i,i,{i),* for Re Q G [Qi^\k), oo) and even k 
for ReQ G [Ql^\k), oo) and odd k 

Therefore, Q = 1, Q = 3 and Q = 5 (only when Q^^\k) > 5^ are isolated limiting points; 
and Q = 2, Q = A, and Q = Q^^\k) are non-isolated limiting points. 

Notice that the above conjecture is also valid for = 3, except that in this case, there 
is no region where the eigenvalue /ifc+i,3,(3),* becomes dominant, since Qf\3) < 4. It IS 
worth noting that the dominant eigenvalue comes from the completely symmetric irreducible 
representation of Si for £=1,2, 3. 

Assuming that the dominant eigenvalue on the real Q-axis invariably comes from the 
completely symmetric representation, it is possible to extend the numerical determination of 
Qc^^fc) to higher values of k. The relevant transfer matrices can be numerically diagonalized 
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k 


Ql^\k) 


1 


3 


2 


3.6180339887 


3 


3.7818423129 


4 


4.5697435537 


5 


4.9029018077 


6 


5.1079785012 


7 


5.2352605291 


Q 
O 


0.oz4DyDDyUo 


9 


5.3886186958 


10 


5.4364766073 


11 


5.4729804532 



Table 2: Values of Qi '{k) for the generalized Petersen graphs G{nk, k) for 1 < k < 11. 



by a standard iterative scheme that uses their decomposition as a product of sparse matrices, 
as in (3.2). The resulting values of Q^f'\k) are shown in Table 2. 

The values of Qc^\k) are seen to be well fitted by a power-law Ansatz: Q^\k) = 
Qc*^ + Ak~^. To prevent possible parity effects, we performed the fits for odd and even 
values of k separately. Both fits were consistent with each other and give a common limit 
Qc*^ = 5.69(1). In fact, we tried a fit of the full data with an Ansatz of the form 

We found that 1) Qi*^ = 5.69(1), 2) that the differences between {A,B) and {A',B') were 
very small, so the parity effects are almost eliminated by fixing a common limit for the odd 
and even data sub-sets, and 3) the fits were rather stable (the fits for data 5 < A; < 9, 
6 < A; < 10, and 7 < /c < 11, gave almost the same results). Finally, we also tried to 
extrapolate our numerical data using the Bulirsch-Stoer algorithm for the odd- and even-Zc 
subsets, and for the full data set. The convergence is rather slow, and our best estimate for 
the limit is ^ = 5.685(10), which agrees with the previous estimates. We conclude that: 

gf) = lim g(*)(A;) = 5.69 ±0.01. (6.5) 

k—^oo 

Formalizing slightly these experimental results, we conjecture that the exists a real number 
^ ^ 5.69 that is an accumulation point of real fiow zeros for the graph family G{nk, k), in 
the limit lim„_j,oo linifc_j.oo) and that no generalized Petersen graph has a real fiow zero above 
g^*^ . The hope that no graph family fares better than the generalized Petersen graph is our 
main motivation — and evidence — for stating Conjecture 1.9 of Ref. [30]: 

Conjecture 6.3 For any bridgeless graph G, ^g{Q) > for Q e [6, oo). 

This conjecture has a great interest in Combinatorics, especially in connection to Tutte's 
five-fiow conjecture [36,55,56]: 
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Conjecture 6.4 For any bridgeless graph G, $g(5) > 0. 



This conjecture seems at least as difficult to prove or disprove as the famous four-color 
conjecture (which eventually turned out to be a theorem [1]). 

7 Chromatic polynomial for the generalized Petersen 
graphs 

We have also symbolically computed the chromatic polynomial PG{nk,k)iQ) for the gen- 
eralized Petersen graphs G{nk, k) with 1 < A; < 6 by specializing the full partition functions 
ZiG{nk,k){Q)V) of Scctious 3 and 4 to the case v = —1. Again, we have checked our symbolic 
results against exact computations made by using maple for the smallest values of n > 1. 

Again, in addition to finding the chromatic roots for finite Petersen graphs G{nk,k), 
we are interested in their accumulation sets as n — )■ oo. We have obtained in this way, for 
each value of k, the corresponding limiting curve of non-isolated limiting points. We 
have used the direct search method [46] to obtain this curves with high precision arithmetic. 
These curves are shown in Figures 11-13. In the cases k = 7,8,9, the chromatic zeros 
were obtained by using maple. The corresponding limiting curves Bk were computed by 
numerically diagonalizing the relevant transfer matrices by the Arnoldi algorithm and the 
use of their decomposition as a product of sparse matrices, as in (3.2). This latter procedure 
was implemented in C. 

Contrary to what happens for the flow-polynomial case, the limiting curves are bounded 
in the complex Q plane (i.e., there are no outward branches going to infinity). This is 
expected because a theorem by Sokal [52, Corollary 5.3 and Proposition 5.4] states that if G 
is a loopless finite undirected graph of maximum degree A, then all the chromatic zeros lie in 
the disk \Q\ < KA, with A ^ 7.963906. For the generalized Petersen graphs A = 3, and the 
bound given by this theorem is far from sharp for this family of graphs. If we restrict to the 
cases with k > 3 (see Figures ll(c)-(d), 12, and 13), we notice some regularities about the 
isolated limiting points and the dominant eigenvalues in the complex Q-plane. If we denote 
by Qc^^ the largest real value of Q where Bk crosses the real Q-axis, then we conjecture that: 

Conjecture 7.1 Fix k >3. Then, the dominant eigenvalue in the regions that intersect the 
real Q-axis is 



Therefore, Q = 1 is the only isolated limiting point; and Q = 0, Q = 2 and Q = Qc (k) are 
non-isolated limiting points. 

Note that, contrary to the fiow-polynomial case, the integer Q = is a (non-isolated) 
limiting point. The similarities of the limiting curves can be appreciated from Figure 14(a), 
where we show together the limiting curves B^ ior 1 < k < 9. The limiting curves seem 
to converge to some infinite- width curve Boo, and the picture is qualitatively very similar 




for ReQ e (-oo,0] 
forReQ e [0,2] 
forReQe[2,QP{k)] 
forReQ e[QPik), oo) 



/^fc+l,2,(2),* 
A*fc+1,0,* 



(7.1) 
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Complex zeros Pq(^^ j) Complex zeros Pg^gn 2) 




-10 12 3 -10 12 3 

Re(Q) Re(Q) 

(c) (d) 



Figure 11: Complex zeros of the chromatic polynomial and limiting curves Bk in the complex 
Q-plane for the generalized Petersen graphs G{nk, k) with k = 1 (a), k = 2 (b), k = 3 (c), 
and k = 4 (d). For each value of k, the zeros correspond to the generalized Petersen 
graphs G{10k, k) (black □) and G{20k, k) (red o). Each region is labeled with the sector the 
dominant eigenvalue belongs to (e.g., i = 1). 
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Complex zeros Pq^^^ 5) Complex zeros Pg^gn 6) 




-10 12 3 -10 12 3 

Re(Q) Re(Q) 

(c) (d) 

Figure 12: Complex zeros of the chromatic polynomial and limiting curves Bk in the complex 
Q-plane for the generalized Petersen graphs G{nk, k) with k = 5 (a), k = 6 (b), k = 7 (c), 
and k = 8 (d). The black squares (□) correspond to the zeros of G(50,5) (a), G(60,6) (b), 
G(21,7) (c), and ^(18, 8) (d). The red circles (o) correspond to the zeros of (j(100,5) (a), 
G(120,6) (b), G(28,7) (c), and G(24, 8) (d). The labels are as in Figure 11. 
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Complex zeros Pq^q^ g) 




-10 12 3 

Re(Q) 



Figure 13: Complex zeros of the chromatic polynomial and limiting curves Bk in the complex 
Q-plane for the generalized Petersen graphs G{9n,9). The black squares (□) correspond to 
the zeros of G{18, 9), and the red circles (o) correspond to the zeros of G{27, 9). The labels 
are as in Figure 11. There are two tiny complex-conjugate oval-like regions (as for the other 
odd values of k) around {Q,v) ~ (2.80,0.36); see Figure 14(b). 



to those for the square-lattice chromatic polynomial with both cyclic [28] and toroidal [29] 
boundary conditions. For odd values of k, the limiting curve displays two complex- 
conjugate regions on its right-most part. Inside these regions, the dominant eigenvalue 
comes from the i = 1 sector. These regions are shown for = 3, 5, 7 in Figure 14(b). It is 
clear that, as k increases, they become smaller and closer to the real Q-axis. Therefore, we 
conjecture that they completely disappear in the limit k — )■ oo. 

Assuming that the dominant eigenvalue on the real Q-axis invariably comes from the 
completely symmetric representation, it is possible to extend the numerical determination 
of Qi^\k) to higher values of k. The resulting values of Qi^\k) are shown in Table 3. 

The values of Qi^\k) show parity effects, so we fitted the data for k even and the data 
for k odd separately. Using the same techniques as for Q^^\ we obtain in this case 

= Mm QP{k) = 2.861 ±0.003. (7.2) 

The power-law fit to the odd-fc data gave a contant term ^ 2.86(1), while for the even-A; data 
the result was ~ 2.86(8). If we used the full data set with the Ansatz (6.4), we obtained 
^ 2.861 for 7 < A; < 11, and ~ 2.864 for 6 < A; < 10. The error bar was taken as the 
difference between the last two estimates. 
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Figure 14: (a) Limiting curves Bk for the chromatic polynomial in the complex Q-plane for 
the generalized Petersen graphs G{nk, k) with k = 1 (black), k = 2 (gray), k = 3 (light blue), 
k = 4 (green), k = 5 (orange), k = 6 (pink), k = 7 (violet), k = 8 (red), and k = 9 (navy 
blue), (b) Zoom for odd values of k = 3,5,7, where we find two small complex-conjugate 
closed regions. Inside them, the dominant eigenvalue comes from the i = 1 sector. Notice 
the change of scale with respect to panel (a). We do not show the curve for = 1, as then 
the small regions for k = 7 would be difficult to see. The labels are as in Figure 11. 



8 Discussion 

In this concluding section we would like to make an unified picture of the results presented 
in this paper, and to propose a number of conjectures that should stimulate further research. 
We have seen that the Q-state Potts model with temperature parameter v comprises both 
the fiow polynomial {v = —Q), and the chromatic polynomial {v = —1) as special cases. 
It is therefore useful to compare our results on the partition function and its fiow- and 
chromatic-polynomial specializations for the particular class of non-planar graphs G{nk, k) 
with what is known for the same objects evaluated on certain classes of planar graphs. Our 
graphs are taken to consist of n identical layers of width L vertices and we impose periodic 
boundary conditions in the n-direction. 

Before starting with the main discussion it is worth recalling the BKW theorem [10-13,53] 
in a simple way: Let D be a domain (connected open set) D in the complex plane, and let 
c^M, /^i, /Um {M > 2) be analytic functions on D, none of which is identically 
zero. Then, for each integer n > 0, define the functions 

M 

fc=i 

Each function has a set of complex zeros, and we are interested in the accumulation set of 
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k 



WW) 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 



2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 



6383423072 
5054257523 
6971127211 
6947536196 
7631358388 
7682556521 
7983688222 
8023368562 
8176338815 
8203525818 



Table 3: Values of Qc {k) for the generalized Petersen graphs G{nk, /c) for 1 < /c < 11. 



these zeros in the limit n oo. Then, under a mild "no-degenerate-dominance" condition, 
these limiting zeros can be characterized as follows: a point z lies in this limiting set if and 
only if either 

(a) There is a unique dominant eigenvalue at z (i.e., > |A*i| for all i ^ k), and 

ak{z) = 0; 

(b) There are two or more dominant eigenvalues /ik at z. 

Let us now consider first the planar chromatic case. The general situation is expected to 
be the following: 

Conjecture 8.1 Consider a family Gn,L of planar connected graphs, consisting of n identical 
layers of width L vertices, with periodic boundary conditions in the n-direction. Let Al (resp. 

Bl) denote the set of limiting points of chromatic zeros for Gn,L, in the limit n — )■ oo, that 
satisfies condition (a) (resp. condition (h)) of the BKW Theorem. Further let Bp he the p-th 
Beraha number (1.12), and let 



B' 



I even 



{Q^Bp-. pe2N} , 



B 



lodd 



p e 2N + 1} 



p 




Then: 



2. For any p ^ N, the dominant eigenvalue of the transfer matrix for 



Q e (i?2p,i^2(p+i))n[0,gf)(L)) 



comes from the sector with p marked clusters. 
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Compelling evidence for Conjecture 8.1 comes from the results of Ref. [28], where it was 
shown that Conjecture 8.1 is correct for cyclic strips of the square and triangular lattices 
with L < 8. The corresponding values of Qi^\L) can be found in [28].^^ 

The natural step suggested by these finite-L results is to go to the infinite-volume limit, 
limL_j,oo- We expect that for each family of planar graphs Gn,L, defined as above, the limit 
Qi^'' — liuiL^^ Qi^\L) exists and satisfies < Qc^^ < 4 [28,49,50]. In particular for cychc 
strip graphs of the square lattice, we have that Qc^^ — 3 [4,32]. For cychc strips of the 
triangular lattice, we expect that either ^ 3.8196 [5,6], or qP = = 2 + [31]. 

The deep meaning of these results is that the special role of the Bcraha numbers in the 
production of real chromatic zeros is brought out whenever the chromatic line (v = —1) 
intersects the BK phase in the {Q,v) plane of the Potts model [49,50]. For each graph 
family, and for each Q e [0,4), it can be argued [33] that the extent of the BK phase 
is a finite non-empty interval v^{Q) < v < v+{Q) with the properties v+{0) — and 
limQ^4 = limQ^4 f In particular, other choices of the temperature curve v{Q) 
are possible and might in some cases lead to higher values of Qc than those quoted above. ^'^ 

The bound on Qc given above is even sharp: i.e., there exists a family of planar graphs 
Gn,L, defined as above, and a curve v — f{Q) in the {Q, v) plane, so that the limiting set of 
zeros for the Potts-model partition function Zc{Q,v) satisfies Conjecture 8.1 and Qc = 4. 
Indeed, this is the case [33,49,50] for cyclic strips of the square lattice with the choice [2] 
v"^ = Q for < Q < 4 and —2 < f < 0. Another example is provided by cyclic strips of the 
triangular lattice with the choice [7] + 3f ^ = Q for < Q < 4 and —2 < v < 0. 

As the amplitudes do not depend on i>, we can extend Conjecture 8.1 to the zeros of the 
full partition function Za(^n,L)iQ,v) in the {Q,v) plane: 

Conjecture 8.2 Consider a family G^^l of planar connected graphs, consisting ofn identical 
layers of width L vertices, with periodic boundary conditions in the n-direction. Let the sub- 
sets of Beraha numbers B^^^^ and B°'^'^ be as in Conjecture 8.1. Then: 

1. There are two curves v±{Q) such that < v^{Q), ^+(0) = 0, and a number 
Qo < 4; such that \imQ^Q^^v+{Q) = limQ_^Q„ v_{Q). 

Furthermore, consider any smooth curve v{Q) in between these two curves: < v{Q) < 

v+{Q) with Q G [0,(5o]- Let us denote its graph in the {Q,v) plane as C. Let Al (resp. Bl) 
denote the set of limiting points of partition-function zeros for G^^l, in the limit n — > oo, 
that satisfies condition (a) (resp. condition (b)) of the BKW Theorem. Then: 

2. For each curve v{Q) with the above properties, there exists, for each L, a number 
Qc{L) e (0, Qo) so that 

(a) Re{AL n C) = B"'^'^ n [0, Qc{L)) 

(b) Re{BL n C) = n [0, Qc{L)) U {Qc{L)} 

^^In Ref. [28] Qi^\L) is simply denoted by Qc{L), as there is no danger of confusion with the similar 
quantity Qc*^ associated to the flow polynomial. 

^^We use Qc in this case to note that the limit n — >^ oo is taken along a curve v{Q). Indeed, when 
v{Q) = -1, Qc = and when v{Q) = -Q, Qc = Qi^^- 
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3. For any p E N, the dominant eigenvalue of the transfer matrix for 

ReQ e {B2p,B2ip+i))n[0,Q,iL)) 

and ImQ in between the curves v±{Q), comes from the sector with p marked clusters. 

Let us finally mention that Conjectures 8.1 and 8.2 enjoys strong support from the special 
representation theory of the quantum group at roots of unity [49,50]. 

Remark. A natural question is to determine the value of the number Qq defined in the 
previoTis conjcctTirc for planar strip graphs. Prom the exact results for the free energy of 
the Q-statc Potts antiferromagnet on the square, triangular, and hexagonal lattices [2-7], 
one would be tempted to claim that Qq = 4. This value is also supported by the numerical 
results for the (4, 8^) and (3, 12^) lattices, obtained using the Jacobsen-ScuUard method [35]. 
However, the same method [34, 35] suggest that Qq < A for the kagome lattice. Therefore, 
we prefer to be conservative and not make any firm conjecture about the value of Qq for 
general planar lattices. 

Turning now to the partition-function zeros of non-planar graphs. Figures 3, 10, and 14, 
Conjectures 6.2 and 7.1, and the observations made in Section 4, makes us confident that 
the statements of Conjecture 8.1 hold true also in the partition-function case for non-planar 
graphs, provided that one replaces even/odd Beraha numbers by even/odd integers. On a 
fundamental level, the replacement of Beraha numbers by integers stands out most clearly 
by comparing the eigenvalue amplitudes in the two cases [see, e.g., [28, Eqs. (2.28)-(2.29)] 
and Eq. (3.4)]. To be precise: 

Conjecture 8.3 Consider a family Gn,L of non-planar connected graphs, consisting of n 
identical layers of width L vertices, with periodic boundary conditions in the n-direction. 
Then: 

1. There are two curves v±{Q) such that V-{Q) < V-^{Q), v+{0) — 0, and there is a value 
Qo (possibly Qq^oo), such that Yhxiq^q^ ^'+(<5) = limg-^Qo 

Furthermore, consider any smooth curve v{Q) in between these two curves: V-{Q) < v{Q) < 
v^{Q) with Q G [0,Qo]- Let us denote its graph in the {Q,v) plane as C. Let Al (resp. Bl) 
denote the set of limiting points of partition-function zeros for Gn,L, in the limit n oo, 
that satisfies condition (a) (resp. condition (b)) of the BKW Theorem. Then: 

2. For each curve v{Q) with the above properties, there exists, for each L, a number 
Qc{L) e (0, Qo) so that 

(a) Re{AL n C) = (2N - 1) n [0, Qc{L)) 

(b) Re(i3LnC) = 2Nn [0,Qc(^))U{(5c(-t')}- For some curves v{Q) the value might 
not be present in this set. 

3. For any p & N, the dominant eigenvalue of the transfer matrix for 

ReQ e (2(p-l),2p)n[0,a(i^)) 

and ImQ in between the curves v±{Q), comes from the fully symmetric irreducible 
representation (p) of the sector with p marked clusters. 
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This conjecture has been vahdated in this paper for Gn,L — G{nk,k), the generahzed 
Petersen graphs, with L = k + 1, and /c < 7 for the flow polynomial v{Q) = —Q, and the 
chromatic polynomial v{Q) = — 1. We have also presented numerical evidence that the limits 
Qc exist along these two curves v{Q). In addition, we have directly checked the consistency 
of this picture (at least for /c < 7) by considering directly the partition-function zeros along 
several lines v{Q) — —pQ, v{Q) — —Q/p, and v{Q) — —p with p a positive integer (see 
Section 4). 

Remcirk. The set Re{BL n C) does not contain the value for the flow polynomial; but it 

contains this value for the chromatic-polynomial case. The difference might be due to the 
fact that when Q < we enter in the ferromagnetic regime in the former case, while we stay 
in the antiferromagnetic or unphysical regime in the latter case. 

To conclude, we recall that for the Potts model on planar graphs Gn,L, when Q — 
(2cos(7r/p)) is a Beraha number the representation theory of the underlying quantum al- 
gebra [39] leads to massive degeneracies in the spectrum of the transfer matrix. In partic- 
ular, depending on p, complete sectors of eigenvalues are contained in other sectors in an 
inclusion-exclusion fashion, completely independent of the value of v. We have seen in Sec- 
tion 5 that a similar phenomenon occurs for the Potts model on non-planar grsiphs when Q is 
a non-negative integer. There are striking parallels between the detailed inclusion-exclusion 
scenarii in the two cases, once again provided that one "translates" between the two cases 
by replacing the Beraha numbers by non-negative integers. 



A Potts— model partition function for the simple-cubic 
graphs 

A natural question is that our conclusions are based on a rather particular family of non- 
planar graphs. In this appendix we want to examine a more "physical" family of graphs: the 
graphs Sc(L, n) formed by a simple cubic graph of section of size LxLxn with cyclic bound- 
ary conditions (i.e., free boundary conditions in the transverse "space-like" two-dimensional 
layers, and periodic boundary conditions in the longitudinal "time-like" direction). The 
chromatic polynomial for this family was previously considered in Ref . [45] . 

The graph Sc(L, n) contains L'^n vertices, 2(L — l)Ln horizontal edges, and L'^n vertical 
edges. If we denote the vertices in Sc(L, n) as points in Z^, the vertex set is 

V{Sc{L, n)) = {{x,y,z)\0<x,y<L-l,0<z<n-l}, (A.l) 

and the edge set is the union of the sets E{Sc{L,n)) — E^U EyU E^, where 

= {{{x,y,z),{x + l,y,z))\0<x<L-2,0<y<L-l,0<z<n-l}{A.2si) 
Ey = {{{x,y,z),{x,y + l,z))\0<x<L-l,0<y<L-2,0<z<n-l}{A.2h) 
E, = {{{x,y,z),ix,y,z + l))\0<x,y <L-1, 0<z <n-l} (A.2c) 

where we have identified z = n and z = 0. 

We can apply the transfer-matrix formalism of Section 3 to the family Sc(L,n). If we 
label the vertices on a horizontal layer as {x,y) with < x,y < L — 1, then the transfer 
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matrix takes the form Tl 



V • H with 



L-lL-2 L-lL-2 



H 




(A.3a) 



y=0 x=0 x=0 y=0 



L-l L-l 



V 



n n 



(A.3b) 



y=0 x=0 



The order of the operators y^(x,y),{x+i,y) and y^(x,y),{x,y+i) in (A. 3a) is of no importance, as all 
the horizontal operators commute. The same also holds for the vertical operators V^^.^^^) in 
(A. 3b). However, the operators H and V do not commute. 

The computation of the partition function for the graph Sc(L,n) is quite demanding: 
the number of partitions we have to deal with is given by the Bell number B2L2 [30, and 
references therein]. This number grows very fast as a function of L. The formula for the 
partition function for G = Sc(L, n) is given in terms of traces of the relevant diagonal bocks 
of the full transfer matrix T^: 



e=o xeSe 

Some of the eigenvalues of the transfer matrices £ ^ may coincide. 

In practice, we have been able to compute only the exact partition function for L — 2. 
The structural properties of the transfer matrix T2 can be summarized as follows: 

• The trivial eigenvalue 112,4 = v'^ appears for all values of the number of links < i < 
— 4 with multiplicities 2, 6, 12, 24. This implies that the corresponding amplitude 



• For i = there are two eigenvalues that appear twice; therefore in the "complete" 
representation of ^sc(2,n), they should have the amplitude 2ao — 2. 



• For i = 1, there are eight eigenvalues that appear twice; therefore their amplitude is 
2ai =2(Q-1). 



• For £ = 2, there are four eigenvalues that appear twice in both representations (2), 
and (1, 1). Its amplitude is a2,{2) + «2,(i,i) = 2((5^ — 3Q + 1). In addition, there are six 
eigenvalues in the representation (2) that appear twice; and another six eigenvalues in 
(1, 1) that appear twice, too. Their corresponding amplitude should be twice the value 
Q;i,a- 

• For £ = 3, the representations (3) and (1, 1, 1) have exactly the same eigenvahics. Their 
amplitude is (Q - 1){Q^ - 5Q + 3)/3. Two of these eigenvalues appear twice in each 
representation, therefore their amplitude is twice the latter value. The representation 
(2, 1) has all its eigenvalues doubled; hence their amplitude is 2Q{Q — 2){Q — 4)/3. 

The number of non-trivial distinct eigenvalues in each £ sector are 12, 27, 39, 14 for £ = 
0,1,2,3, respectively. We then obtain a "complete" decomposition of Zsc(2,n) in terms of 
distinct eigenvalues Hl,(.,\,s- The full expression is rather cumbersome, so we refrain from 
writing it explicitly here. 

^^Please, do not confuse this Bell number with a Beraha number (1.12). 




(A.4) 



is 74 = Q2 _ gg3 ^ 2og2 - isg + 1. 
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Real zeros of Z, 



'Sc(2,n) 



o 




1=0 



> 







2 



4 



6 



Q 



Figure 15: Real zeros of the Potts-model partition function and limiting curves B2 in the 
plane {Q,v) for the simple-cubic graphs Sc(2,?t,) for n = 9 (black □), and n = 10 (red o). 
The rest of the notation is as in Figure 1. 



We have computed the zeros of the partition function Zsc{2,n){Q, v) along the same lines 
in the real {Q, v) plane as for the generalized Petersen graphs. Figure 15 shows these zeros 
for the simple-cubic graphs Sc(2,9) and Sc(2, 10). 

We have also computed the limiting curve B2 in the real {Q , v)-p\a.iae; this is depicted 
in Figure 15. The phase structure for this family of graphs is similar to that described for 
the generalized Petersen graphs. In particular, we can observe two phases enclosed by the 
limiting curve B2 and characterized by £ = 1 and i = 2 respectively. The former contains 
the interval Q G (0,2), and the latter, the interval Q G (2,4). In each phase, the dominant 
eigenvalue comes from the fully symmetric representation of the corresponding group Sg. 
Therefore, as the expressions for these eigenvalues are the same as for the generalized Petersen 
graphs, the values Q = 1 and Q = 3 are isolated limiting curves. Notice that there is no 
isolated limiting point at Q = -B5. 

Finally, we observe close to the lower boundary of the BK phase, two regions with domi- 
nant pairs of complex conjugate eigenvalues. This feature was absent in the Petersen graphs. 

In conclusion, we find a phase diagram for this family of non-planar graphs which is 
qualitatively similar to that of the Petersen graphs. This supports our belief that the results 
we have found for the Petersen graphs are general for at least all families of recursive non- 
planar graphs. 

Acknowledgments 

We thank Gordon Royle for suggesting us to study the generalized Petersen graphs, and 
Alan Sokal for many illuminating discussions. 



43 



Both authors also thank the Isaac Newton Institute for Mathematical Sciences, Univer- 
sity of Cambridge, for hospitality during the programme on Combinatorics and Statistical 
Mechanics (January-June 2008), where this project started. J.S. also warmly thanks the 
Ecole Normale Superieure for hospitality in June 2009, June 2011, and December 2012. 

The research of J.L.J, was supported in part by the Agence Nationale de la Recherche 
(grant ANR-lO-BLAN-0414: DIME) and the Institut Universitaire de Prance. The research 
of J.S. was supported in part by Spanish MINECO grants FIS2012-34379 and MTM2011- 
24097, and by U.S. National Science Foundation grant PHY-0424082. 

References 

[1] K.l. Appel and W. Hakcn, llhnois J. Math. 21, 429-567 (1977). 
[2] R.J. Baxter, J. Phys. C: Solid State Phys. 6, L445-L448 (1973). 

[3] R.J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, London- 
New York, 1982). 

[4] R.J. Baxter, Proc. Roy. Soc. London A 383, 43-54 (1982). 

[5] R.J. Baxter, J. Phys. A: Math. Gen. 19, 2821-2839 (1986). 

[6] R.J. Baxter, J. Phys. A: Math. Gen. 20, 5241-5261 (1987). 

[7] R.J. Baxter, H.N.V. Temperley and S.E. Ashley, Proc. Roy. Soc. London A 358, 535-559 
(1978). 

[8] A. Bedini and J.L. Jacobsen, J. Phys. A: Math. Theor. 43, 385001 (2010), 
arXiv:1003.4847. 

[9] S. Beraha, Infinite non-trivial families of maps and chromials. PhD Thesis, John Hop- 
kins University, 1975. 

[10] S. Beraha and J. Kahane, J. Combin. Theory Ser. B 27, 1-12 (1979). 

[11] S. Beraha, J. Kahane and N.J. Weiss, Proc. Nat. Acad. Sci. USA 72, 4209 (1975). 

[12] S. Beraha, J. Kahane and N.J. Weiss, Limits of zeroes of recursively defined families of 
polynomials, in Studies in Foundations and Combinatorics (Advances in Mathematics 
Supplementary Studies, Vol. 1), ed. G.-C. Rota, pp. 213-232 (Academic Press, New 
York, 1978). 

[13] S. Beraha, J. Kahane and N.J. Weiss, J. Combin. Theory Ser. B 28, 52-65 (1980). 

[14] D.A. Bini and G. Fiorentino, Numerical computation of polynomial roots using MPSolve 
version 2.2 (January 2000). Software package and documentation available for download 
at f tp: //f tp.dm.unipi . it/pub/mpsolve/MPSolve-2. 2 . tgz 

[15] D.A. Bini and G. Fiorentino, Num. Algorithms 23, 127-173 (2000). 



44 



[16] G.D. Birkhoff and D.C. Lewis, Trans. Amer. Math. Soc. 60, 355-451 (1946). 

[17] S.-C. Chang, J.L. Jacobsen, J. Salas and R. Shrock, J. Stat. Phys. 114, 763-823 (2004), 
arXiv:cond-mat/0211623. 

[18] S.-C. Chang, J. Salas and R. Shrock, J. Stat. Phys. 107, 1207-1253 (2002), arXivxond- 
mat/0108144. 

[19] S.-C. Chang and R. Shrock, Physica A296, 131-182 (2001), arXiv:cond-mat/0005232. 

[20] S.-C. Chang and R. Shrock, J. Stat. Phys. 112 815-879 (2003), arXiv:cond- 
mat/0205424. 

[21] S. Dasmahapatra and P. Martin, J. Phys. A: Math. Gen. 29, 263 (1996). 

[22] CM. Fortuin and P.W. Kasteleyn, Physica 57, 536-564 (1972). 

[23] A. Guttmann and I.G. Enting, J. Phys. A 27, 5801-5812 (1994), arXiv:hep-lat/9312083. 
[24] T. Halverson and A. Ram, Eur. J. Comb. 26, 869-921 (2005). 

[25] B. Jackson, J. Geom. 76, 95-109 (2003), arXiv:math/0205047. 

[26] J.L. Jacobsen, J.-F. Richard, and J. Salas, Nucl. Phys. B 743, 153-206 (2006), 
arXi v: cond-mat /0511059. 

[27] J.L. Jacobsen and J. Salas, J. Stat. Phys. 104, 701-723 (2001), arXiv:cond- 
mat/0011456. 

[28] J.L. Jacobsen and J. Salas, J. Stat. Phys. 122, 705-760 (2006), arXiv:cond- 
mat/0407444. 

[29] J.L. Jacobsen and J. Salas, Nucl. Phys. B 783, 238-296 (2007), arXiv:cond- 
mat/0703228. 

[30] J.L. Jacobsen and J. Salas, Is the five-flow conjecture almost false?, arXiv:1009.4062. 

[31] J.L. Jacobsen, J. Salas and A.D. Sokal, J. Stat. Phys. 112, 921-1017 (2003), arXiv:cond- 
mat/0204587. 

[32] J.L. Jacobsen and H. Saleur, Nucl. Phys. B 743, 207-248 (2006), arXivxond- 
mat/0512058. 

[33] J.L. Jacobsen and H. Saleur, J. Stat. Phys. 132, 707-719 (2008), arXiv:0803.2665. 

[34] J.L. Jacobsen and C.R. ScuUard, J. Phys. A: Math. Theor. 45, 494003 (2012); 
arXiv:1205.0622. 

[35] J.L. Jacobsen and C.R. ScuUard, J. Phys. A: Math. Theor. 46, 075001 (2013); 
arXiv:1211.4335. 

[36] F. Jaeger, Nowhere-zero flow problems, in L.W. Beineke and R.J. Wilson (eds.). Selected 
topics in graph theory 3 (Academic Press, 1988). 

45 



[37] W. Janke and R. Villanova, Nucl. Phys. B 489, 679-696 (1997), arXiv:hep-lat/9612008. 

[38] J. Oxley and D. Welsh, Combin. Probab. Comput. 11, 403 (2002). 

[39] V. Pasquier and H. Saleur, Nucl. Phys. B 330, 523-556 (1990). 

[40] V. Pasquier and H. Saleur, Nucl. Phys. B 330, 523-556 (1990). 

[41] R.B. Potts, Proc. Camb. Phil. Soc. 48, 106-109 (1952). 

[42] N. Read and H. Saleur, Nucl. Phys. B 613, 409 (2001), arXiv:hep-th/0106124. 

[43] J.-F. Richard and J.L. Jacobsen, Nucl. Phys. B 750, 250-264 (2006), arXiv:math- 
ph/0605016. 

[44] J.-F. Richard and J.L. Jacobsen, Nucl. Phys. B 769, 256-274 (2007), arXiv:math- 
ph/0608055. 

[45] J. Salas and R. Shrock, Phys. Rev. E 64, 0111111 (2001), arXiv:cond-mat/0102190. 

[46] J. Salas and A.D. Sokal, J. Stat. Phys. 104, 609-699 (2001), arXiv:cond-mat/0004330. 

[47] J. Salas and A.D. Sokal, J. Stat. Phys. 135, 279-373 (2009), arXiv:0711.1738. 

[48] J. Salas and A.D. Sokal, J. Stat. Phys. 144, 1028-1122 (2011), arXiv:1002.3761. 

[49] H. Saleur, Comm. Math. Phys. 132, 657-679 (1990). 

[50] H. Saleur, Nucl. Phys. B 360, 219-263 (1991). 

[51] R. Shrock, Discrete Math. 231, 421-446 (2001), arXiv:cond-mat/9908387. 

[52] A.D. Sokal, Comb. Prob. Comp. 10, 41-77 (2001), arXiv:cond-mat/9904146. 

[53] A.D. Sokal, Comb. Prob. Comp. 13, 221-261 (2004), arXiv:cond-mat/0012369. 

[54] W.T. Tutte, Proc. London Math. Soc. 51, 474-483 (1950). 

[55] W.T. Tutte, Can. J. Math. 6, 81-90 (1954). 

[56] W.T. Tutte, Graph theory (Addison- Wesley New York, 1984). 

[57] M.E. Watkins, J. Combin. Theory 6, 152-164 (1969). 

[58] F.Y. Wu, Rev. Mod. Phys. 54, 235-268 (1982) 235; errata Rev. Mod. Phys. 55, 315 
(1983). 

[59] F.Y. Wu and Y.K. Wang, J. Math. Phys. 17, 439 (1976). 

[60] C.-Q. Zhang, Integer flows and cycle covers of graphs (Marcel Dekker, New York, 1997). 



46 



